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Abstract: In this paper, a new algorithm based on differential geometry view- 
point to solve the 3D rotating Navier-Stokes equations with complex Boundary is 
proposed, which is called Bi-parallcl algorithm. For example, it can be applied to 
passage flow between two blades in impeller and circulation flow through aircrafts 
with complex geometric shape of boundary. Assume that a domain in can be 
decomposed into a series sub-domain, which is called "flow layer" , by a scries smooth 
surface Qk,k = 1,...,M. Applying differential geometry method, the 3D Navier- 
Stokes operator can be split into two kind of operator: the "Membrane Operator" on 
the tangent space at the surface 3^ and the "Bending Operator" along the transverse 
direction. The Bending Operators are approximated by the finite different quotients 
and restricted the 3D Naver-Stokes equations on the interface surface SS^, a Bi-Parallel 
algorithm can be constructed along two directors: "Bending" direction and "Mem- 
brane" directors. The advantages of the method are that: (1) it can improve the 
accuracy of approximate solution caused of irregular mesh nearly the complex bound- 
ary; (2) it can overcome the numerically effects of boundary layer, whic is a good 
boundary layer numerical method; (3) it is sufficiency to solve a two dimensional 
sub-problem without solving 3D sub-problem. 
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tating Navier-Stokes equations; Bi-Parallel Algorithm. 
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1 Introduction 

As well known, numerical simulation for 3D viscous flow in turbo-machinery and circulation flow 
through aircraft meet three kinds of difficulty cause of nonlinearity, three dimensional grid and 
boundary layer effects with complex boundary. In order to overcome last two difficulties, a new 
dimensional slitting method based on differential geometry method is proposed. 

As well known that classical domain decomposition method is that the 3D domain is made in the 
sum of severest overlap or non overlap subdomain, an approximated solutions can be established 
then by solving 3D subproblcm in each subdomain, sec e.g., [2,4,19,24] . The method proposals 
by authors here , called "bi-parallel algorithm ", for 3D Navier-Stokes equations, is that the 3D 
domain fl occupied by the fluids in Jf^ is decomposed into the sum of servery sub domains ( called 
layer) ^ by servery 2D surfaces(2D manifold) 9^, i = 1, 2, • • • , m. 3D Navier-Stokes operators in 
the layer ni_-j^L>fll'^^ under a new coordinate based on the surface Qi can be represented as the sum 
of "membrane" operators on tangent space at and normal (bending) operator to 9. By Eular 
central difference quotient instead of bending operator and restricting 3D Navier-Stokes equations 
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on the Si, a three eomponents and two dimensional Navier-Stokes equations (is called 3C-2D 
NSE) on are obtained. After successively iterations, an approximate solution of 3D NSE can 
be established. It is obvious that the method is different from the classical domain decomposition 
method because we only solve a two dimensional problem in each sub domain without solving 3D 
subproblems. In addition, other advantages of this method are the foUowings: 

(i) For the complex boundary geometry, for examples, in turbo-machinery flow with complex 
shape of blades of impeller, in geophysical flow with real surface of the earth and in the circulation 
flow passing through complex aircraft etc., 80 percent of the freedoms of 3D mesh should be 
concentrated on a thin domain of the boundary layer, even using different methods or finite element 
method, there exists a lot of difficulty and even influence of the accuracy ; in our method the 
distance between two surface can be very small as you wanted, owing to parametrization of surface, 
the domain for 2D-3C subproblem is a bounded domain in 5R^, and the mesh is 2D-mesh; 

(ii) Since interface surfaces arc chosen such that most of flux flows along the interface while 
small flux of the fluid run cross the interface only, iterative method is very well suitable for the 
physical properties of flow; 

(iii) This is a better method for treatment of boundary layer phenomenons, which can refer to 
[19] and references therein. Indeed this is a good boundary layer model and associated algorithm. 

This paper is organized as follows: in section 2 we give the preliminary for the geometry of 
the blade's surface and constriict a new coordinate system; in section 3, we derive the rotating 
Navier-Stokes equations in the new coordinate system; in section 4, we study The equations for the 
average velocity along the rotating direction; in section 5, provided the equations of the Goteaux 
Derivatives of the solutions of Navicr-stokes equations with respect to shape of boundary; in section 
6 we provide 2D-3C Navier-Stokes equations on the Surface 9j ; in section 7, we provide a corrected 
equation for the pressure in order to improvement of the accuracy of the pressure's computation; in 
section 8, we present the Bi-parallel Algorithm; in section 9, we prove the existence of the solutions 
for 2D-3C Navier-Stokes equations and the solution of correcting equation for the pressure, and 
discuss the dependence of the corrected pressure upon the velocity; and finally, an appendix is 
supplied. 

2 Preliminary— The Geometry of the Blade's Surface 

we use Einstein summation convention throughout this paper, in order to distinguish, Greek indices 
and exponents belong to the set {1, 2}, Latin indices and exponents belong to the set {1, 2, 3}, and 
the repeated index implies that we are summing over all of its possible values. The dot product 
and the cross product of two vectors a,b G are denoted hy a ■ b , a x b, respectively, and the 
Euclidean norm |a| = y/a ■ a. 

Let D be an open bounded connect subset of B?, whose boundary 7 is Lipschitz-continuous. 
Suppose D is locally on one side of 7. In the following of this paper, we also suppose that the blade 
of impeller is so thin that it can be described as an smoothing injective mapping 3 = 3?(£') in 
mathematics. In this case, D is the projection area of blade 9 on the meridian plane of impeller. 
Let X = {x") is an arbitrary point in the set D, and da = d/dx" . If the mapping ^{x) is smoothing 
enough, for example, ^{x) G C^(D;B?), then there exist two vectors ea{x) = daJR-ix), which are 
linear independent at all point x G D. ea{x) can be regards as the basis vectors of the tangential 
surface at this point. Meanwhile, the unit normal vector at the same point is 

ei X 62 

n = r. 

|ei X 62! 

So {ea,n) constitute the covariant basis vectors at the point We can further define contra- 

variant basis vectors (e" , n) which satisfy the following relations 

e"-ej}=5'p, e" • n = 0, n n = l, 
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(a) Impeller (b) Flow Passage 

Figure 1: Impeller and Flow Passage 



where Sp denotes the Kronecker delta. It is well known that e" are also on the tangential plane 
to 3 at point SR(a;). 

The covariant components a^^ and contra-variant components a"^ of the metric tensor of 3, 

the Christoffel symbols F^^ , and the covariant components bap and mixed components of the 
curvature tensor of 3 are defined as, respectively, 

* 

bap ■■= n ■ dpBa, bll = a^'^baa- 

The area element along 3 is ^adx, where a = det(aQ,^), ^/a = |ei x e2|. 

Nest, we assume that the impeller is rotating along z-axis with angular velocity uj — (0, 0,w). 
Let {er,eg,k) be the cylindrical basis vectors established on the impeller and rotating with the 
impeller(see Appendix). Let constant N be the number of blade and e = tt/N. Then by rotating 
^ degrees, one blade is rotated to the location of the next one. The flow passage il^ of the impeller 
is the inner part of the boundary dfl^ = Fi„ U Tout U Ft U F^ U 3+ U Further, the blade is 
3 = 3fJ(-D), and an arbitrary point ^(x) on the blade 5 can be expressed as 

m:{x)^x^er + x^e{x)eg + x^k, (2.1) 

where x — (2:^,2;^) £ D are Gauss coordinate system on 5, and Q E C^{D,R) is an smoothing 
enough function. 

It is easy to prove that there exists a family of surfaces 5^, only depending on parameter ^, 
which cover the domain fl^ by mapping 5R(a;; ^) : D — )■ 3^, where 

sR(a;; ^) = x^Br + x^{e^ + e{x))ee + x^k. (2.2) 

It is easy to prove that the metric tensor aap of is homogenous, nonsingular and independent 
of ^, which is given as follows 

aap = ^^^^^^=Sap+r^Mp, a-^ap, = S^, (2 3) 

a = det(a„^) = 1 + r2(e? + ei) > 0. 

Let {x^ , , x'^ ) = (r, 0, z) , the corresponding metric tensor in are {gi'v — 1,52'2' = 
^^j53'3' = ^^9i'j' = 0^*' 7^ j')- According to rule of tensor transformation under coordinate 
transformation we have the following calculation formulae 

dx^ dx^ 

5U = g^r g^j ■ 
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Substituting (2.3) into the above formula, the covariant and contra-variant components of the 
metric tensor of in the new curvihnear coordinate system are give by 



(2.5) 



/ 9a0 = dap-, g?,B = 3/33 = er'^Qp, 933 = ^^r"^ , 9 = det(5ij) = e^r^, 

where [VOp = 9? + Ql and e„ = 

Tensor calculations show that the following proposition is right (see Appendix ), 

Proposition 2.1. let {ea,n) denote the basis vectors in new coordinate system (x,^) , and a 
vector V e can be wrote as v = v°'ea+v'^n. Further more, we have the following representation 
formula in new coordinate system, 
1. Angular velocity vector u: 



1 2 
W = OJ, UJ 



0, uj^ = -we-^Gi, 



(2.6) 



2.Coriolis Force 

■ C = 2u}xw = C'Bi 

< 



U.2 



Ci=0, = -2LurU{w,e), C3 = 2we-i(re2n(t(;,e) + — ), (2.7) 

r 

3. Unit normal vector to 

n= p-^^-^ =n'e„ = -rSx/y/a, = {ery^^/a, (2.8) 
|ei X 62! 

4- Second fundamental form (curvature tensors for 2D manifold) 

bn = 7s(Q2(aii - 1) + rOn), &12 = ^(^1(^12 + ^612), 
^22 ^(62(022 + 1) + re22), b = det{bap) = &11&22 - bj^- 

5. Mean Curvature H and Gaussian Curvature K 

^ = 2^(®2(a + aiia22 + On - 022) + r(a220ii + 011622 - 2012612)), 
K=l = det{b^0)/a. 

It is obvious that each ^ = const corresponds to a surface S>j which has the same geometry 
properties with 3. It is well known that the geometry properties of is completely determined by 
(oa/s) and {bajj) in the following meaning. 

Let denote the set of all third-order orthogonal matrices, and that 0\ = {Q &0^, det((5) = 
1} denotes the subset of all proper third-order orthogonal matrices. Then J+(a;) = C + Qo-x.\sa, 
proper isometry of , where C G R^,Q G O^. We have 

Theorem 2.1. Two immersions G C^{D;R^) and^2 G C^{D;R^) share the same fundamental 
forms (aa/3), (ba^) over an open connected subset D G R^ if and only if 

3*2 = J+o3fii, (2.11) 

Furthermore, If two matrices fields (a^^) e C^(i?;<S^), (bap) € C^(£>;<S^) satisfy Gauss-Codazzi 
equations in D, i.e.. 



(2.9) 



(2.10) 



dpbaa - d^bafS + T^^fe^^ - T^f^b^f, = 0, 

where ^ap,T = ^{^aaar + daa^T — dT-aajs), = a'^'^Taff^T ■ Then there exists an immersion 
3?e C^{D]R^) such that 



\diU X 923?! 
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Because the surface is obtained by a degree rotation of 3, so by using theorem |2.1[ 
we have the surface 9^ have the same geometry properties with 5 for any ^ G [—1, 1], i.e., their 
corresponding geometric quantities Oq/j, bajs, K, H, ■ ■ ■ are same. 

In subsequent discussion, we will often employ the third fundamental tensor of 3, i.e., 

Cap = a^^ba^xbpa, (2.12) 
and inverse matrix (c"'^) — (cajp)^^, (b"^) — (bap)^^ satisfy the following relations, 

b'^^bpx^d^, 2°%A = <5?. (2.13) 
Furthermore, let introduce permutation tensors in Euclid space 

(i,j,k) is even permutation of (1,2,3), 
Eijk = { -y/g, eijk = { (i,j,k) is odd permutation of (1,2,3), (2.14) 

0, otherwise, 

where g = det(g.y). Similar permutation tensors on 2D manifold 9 

(a, P) is even permutation of (1,2), 
£ap — { —y/o-, — { (o^i/^) is odd permutation of (1,2), (2.15) 

0, otherwise. 

3 Rotating Navier-Stokes Equations in the New Coordinate 
System 

we consider the rotating impeller with rotating angular velocity a; ~ (0,0,a;). Under the rotat- 
ing cylindrical coordinate established on the impeller. The motion of fluid in the flow passage is 
governed by the three-dimensional rotating Navier-Stokes equations, i.e., 

dp 

^ + div(pu;) = 0, 
pa = div(T + f, 

dT ' 
P^vi~Q^ + w^^jT) - div(Kgradr) + pd\vw - $ = /i, 

p^p{p,T), 

where p the density of the fluid, w the velocity of the fluid, h the heat source, T the temperature, k 
the coefficient of heat conductivity, Cy specific heat at constant volume, and p viscosity. Further- 
more, the strain rate tensor, stress tensor, dissipative function and viscous tensor are respectively 
given by 

ey(«^) = ^(V.u., -I- V,u.,), e'^iw) = g'^g^^CkUw) = ^{V'w^ + V^w'), 

a'^{w,p) = A^^^-^Ckmiw) - g'^p, A'^''"' = Xg'^g''"^ + p{g'''g^"' + g^'^g^^) (^'^^ 
$ = A^J^-™e,,(«;)e,j(-u;), A = -|/^, 

where gij , and g^^ are the covariant and contra- variant components of the metric tensor of three- 
dimensional Euclidian space in the curvilinear coordinate (a;,^) define by (2.4), respectively. The 
covariant derivatives of velocity vector and Christoffel symbols are 
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The absolute acceleration of the fluid is given by 

dw 



+ {w\')w + 2u} X w + iv X {u} X TZ), 



dt . (3.4) 



where IZ is the radium vector of the fluid particle. The flow passage occupied by the fluids is 
denoted by fJg, which boundary 9f2e is the union of inflow boundary r.i„, outflow boundary Tout, 
positive blade's surface , negative blade's surface 3_ and top wall Ft and bottom wall Ft,, i.e., 

9rje-F = r„uro„tU3_u3+uFturb (Fig.2) (3.5). 

The boundary conditions are given by 



(3.6) 



cr*J(u;,p)nj|r,„ = gj„, (T'J(«',p)nj|r„„, = glut, 
^ + A(T - To) - 0, A > 0. 

We also supply the initial condition 

If the fluid is incompressible and flow is stationary, then the governing equations are 
divii; = 0, 

{wV)w + 2u) X w + Vp- vAYv{e{w)) = -{ujfTl + /, 

w\t„=Q, (3.7) 
(-pn + 2i/e(u;))|r„ = g,;„, 
. (-pn + 2j/e(«;))|r„„t =£/o«t, 

where Fq = 9+ U G_ U Ft U F^. For the polytropic ideal gas and flow is stationary, system (3.1) 
can be wrote as the conservation form 

div(pi(;) = 0, 

dW{pw ®w) + 2puj XW + RS/{pT) = pAw + (A + //)Vdivu; - piuj^TZ, (3.8) 
div[p(i^ + CyT + RT)w] = kAT + XdW{wdivw) + p,div[wWw] + f Ajiop, 

while for isentropic ideas gases, it turns 

div{pw) — 0, 



div(/9t(j <g)w) + 2poj X w + a\7{p^) = 2/idiv(e) + AVdivio - p{ujy7l 



(3.9) 



where 7 > 1 is the specific heat radio and a a positive constant. 

Furthermore, we give the expressions of the power /(3, ^(^y)) done by the impeller and global 
dissipative energy J(9,u;(9)), respectively 

/(3,w(3))= / [ a-n-egujrdQ, J(5,u;(3))= / / / ^MdV, (3.10) 

Under the new coordinate system (2.2), from the discussion in section [2j we know there exists 
mapping between the fixed domain fl — D x [—1,1] and the flow passage Sl^. In the subsequent 
paragraph, we suppose D E is surrounded by four arcs AB, CD, CB, DA, and 

5£) = 7oU7i, 7o = ABUCB, ji^CBUDA. 
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(a) Blade (b) Projection area 

Figure 2: Blade and Projection area D of blade on the meridian. 



There exist four positive functions 70(2;), 70(2), 71 (z), 71 (z) such that 
r := a;2 = 7o(a;i) = 7o(z) q^^^^ ^2 



r := X 



7i(a;i) = 7i(z) o^DA, 



iCD 



fo < 70(2) < ri on 7-0 < 7o(z) < ri onCD, 
fo < 71 (-2) < J"! on DA, tq < 71(2:) < ri oniJC'. 

The corresponding three-dimensional region is expressed as 



dn = ToU Ti,^ Ti = ro„t u r„„ To = r,,^u u - i} u = -i}, 
<9i? = 7oU7i, 7o = (i^nffc)upnft), 71 = puf„„t)u(/?uf„), 

where 3fJ(-) is defined by (2.1) 

We introduce the following Sobolev space, 

Vin) = {v\v e H\n)\ v\^^^ = 0}, H^ifl) = {q\,q € H\n), = 0}. 



(3.11) 



(3.12) 
(3.13) 

(3.14) 



equipped with usual Sobolev norm 
problem (3.7) and (3.9) 



Consider the variational formulation for Navier-Stokes 



Find to e V{n),pe ^^(f^), such that 

a{w, v) + 2{u: X w,v) + b{w,w,v) ~ {p,divv) ^< F,v >, WveV{Q), (3.15) 
{q,divw) = 0, Vg G L'^in), 



and 



Findu; £ V{n),p£ ^^(ri), such that 

a{w, v) + 2(w X w,v) + b{pw, w, v) + (— p + Adivto, divt>) =< F,v >, Vv G V{fl), 
{Vq,pw))^<pw-n,q>\r,, yqeH^iQ), 



where 



<F,v >=< f,v> + <g,v>f^, <g,v >=< gri„, t> 



> ? 



< 9out,V > If 



ro„t 



(3.16) 



(3.17) 
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In order to rewrite equations (3.7) and (3.9) in the new coordinate system, we have to consider 
covariant derivatives of the vector field. Therefore, we first give out the second kind of Christoffel 
symbols in the new coordinate system in terms of 9(see Appendix), i.e.. 




7' 



pa 

^ 3/3 



r^3 = r-i52a + re2e„ r^3 = -£V52a, r33 = £re2, 



(3.18) 



then covariant derivatives of the velocity field is VjW^ = + T-lf.w^, which can be specific 

expressed in the next lcmma(sec Appendix), 

Lemma 3.1. Under the new coordinate system (x^, x^, ^) defined by (2.4), the covariant derivatives 
of the velocity field can be expressed as 



< 



V„«;=^ = e-^x^yw^e^ + e-'w^e^0 + {ex^)-''a2a^{w, 6), 

Vau;" = ^ - x^e62cll{w, 9), V^w^ = + ^ + x^Q2U{w, 6), 



divw=^ + 



+ 



u{w,e) = sw^ + wpef} 



A simple calculation show that the strain tensor can be rewrite in the splitting form 

eij{w) = (t)ij{w) +il)ij{w,Q), 
where the first terms is independent of 6, that is 



(3.19) 



(3.20) 



1 ,dw 



l,dw° 



M = ^(^ + ^). '/'3a(«^) = ^(^+eV 



2 ' dxP 



2^ 



dx°' 



), </)33(t«)=£V( 



-)• (3.21) 



While the second terms contains 6, ^ij{w, 6) = ipij{w)Qx + tl)]^-' {w)Q \@ „ + 6), where 



1 2(9w^^ 

2^'' ^dx^^^^ ' dxf^' 

1 ^,dw ^\,dw^ 2 , \ / N 

1 2r9w^ - dw^ - 2 2 c- c X 



, di 



er 



(3.22) 



2 n/^ ^aai 



^33"(t«) = 0. 



and 



ir2w-5,(e„e^), eL(«;) 



ea/3(«^,0) 
The proof is omitted. 

The following notations are frequently used in the later, 



ler'^w'^Qaa, 



A = 



+ 



divw = 



eUw) = 0. 



(3.23) 



dx" 



For the sake of simplicity, we just consider incompressible flow. Taking into account (3.18), (3.19), 
in the new coordinate system the Navier-Stokes equations can be written in the form, 

Theorem 3.1. Suppose that the blade surface is smooth enough, that is Q is smooth enough, for 
example, 6 S C^{D), then the rotating Navier-Stokes equations in the new coordinate are given 
by _ 

9^ + 9^ + ^= 1 + 9^ = dw,W ^ 



r dx^ 



0, 



J\f''{w,p,e) := -vKw^ - v(re)-'^a^ - i^P^'^S)'^ 

-vPf{Q)S^ - vq'^{Q)w^ +g^P{@)Vfip + g''\&) 
+C''{w,u}) + N^{w,w) = f^, VA; = 1,2,3, 



(3.24) 
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where C{w,oj) is Coriolis forces defined in (A. 1. 4), and 



„3„„3 



rS2^uiw,e)u{w,e), 



+{re)-^U{w, e)(2u>2 + r'^&2^{w, 9)), 



(3.25) 



N''{w,w) 
B''{w,w) 

(9) 
Pf(9) 

ai(0) 



+ d^iw^w^) + n^jw'w^ = + B\w, w), 

lS02, Ps^iO) = -2reS2xe0, 

(r£)-i(<52^9„ + 2r9„^), P^^ = Hp2, 

-[{re)-\Sa,xe2 + 2<52a9A) +£-i5aAA9], P3"3 = -2r-^S2a, 

2e-2(r-352.-9;39^,), 

£-i(r92|V9|2- AO), 

= 2(52„[(5ja|V9|2 - 0629^ + r9^A9 - rexSx^] - r-^52j2oc, 
= (52o(rA9-2a92)e, 

= (r£)-i[r-i(l + a(a22 - 1))9^ + 292<.] + e-^d^AQ, 
- ae2&2, 



(3.26) 



9 = ^ 9 



A9 = 9a„ = 9ii + 922, |V9|2 = 92 + 0i, 



and 



Kai^) = -r5..2&Xea + r-^Sx2Saa, 



^A3(e) - '^fAl©) = -reS2o.ex, TT^s = -r£252„, 

7rL(9) = £-19a. + {re)-'ex{S2a + 02.), 
^ki®) = ^lai^) = ^"'«2A + r-^52x, Trfg = r£92. 

Proof:The Proof see Appendix. 

Let introduce the inner product in the Sobolev space V{Q,) or V{D) 

{w,v) = JlgijW^v^ y/gdx.d£, = /[aa^w"v^ +r^£9/3(w^v^ +w''v^) + r^£^w^v^]r£dxd^, 

a Q 
{w,v)d = Jiaa^w^v^ + r'^sQi3{'w^v^ + w^v^) + r'^e'^w^v^]redx, 

D 

The subscript "D" will be omitted if there is no misunderstanding. 

Next we consider the variational formulation for (3.24) in the new coordinate system, 
into account (2.5), let set 

A{w, v) = {gijN\w,p, 9), v^) = iga^Af" + sr^e^Af^, t;^) + {sr^QpAf^ + e^r^Af^, v^) 
By using index reduction, lift and descent of tensor, we get 

gmkg'"^ = ^l^n^ grnkg''^ = <5m, Cm{w, U}) = gmkC" {w , Lj) , 

PLji^) = 9mkPjK^)y 9mj(9) = gmkqj{Q), Bmiw,w) = gmkB''{w,w). 

Let adopt the notations 

Em{w) = gmk^w^ = dxigmkdxw^) - dxgmkdxw'', 
Nm{w,w) = g,„aN°'{w,w) + g^3N^{w,w) 

= 9mk C'%l' ^ +dfs{w''w^) +TT^jW'w^) = g^k^^^ + Bra{w,w), 
Bm{w,w) = gmk{d^iw''w^) + 'K^jW'w^), 



(3.27) 
(3.28) 

(3.29) 

□ 

(3.30) 

Taking 
(3.31) 

(3.32) 



(3.33) 
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Then 

-,.P,% (6) - uqr., {Q)w^ + SiVpp + Sld^p (3.34) 
+C„i(«;,u;) + Bm{w,w) + Qmkddw'^w'') = fm, 

Remark 3.1. Obviously we have 

{C{w,U!),w) = Cp{w,uj)w^ + C3{w,uj)w^ = 2ru!{uPQ0 - S2j3ll{w,Q))'W^ + 2rsojw'^'w^ 
= 2ru}{w'^w^ldj3 - w'^{ew^ + w^Qx)) + 2reujw^w^ = 0, 

(3.35) 

which is coincide with 2u} x w ■ w = 0. 
Since 

= -dx{vgmkv"^d\w^) + i^g,nkdxw''dxv"^ + vdxgmkdxw^v"^ , 
J Ji-iygmk^w'^v"' + 6^di3pv"']d^d^ 

= / li-dxiiygmkv^dxw'^) + d^{v^p) + vgmkdxw^dxv"^ + i^dxgmkdxw'^v"' - pd0V0]d^d^ 

= / I [-i^9mkv"'dxw''nx+pnpv'^]dsd^ 
{=-iaD 

+ / J[i^gmkdxw''dxv"' + i'dxgmkdxw''v"' -pd/3vi^]dx.d^, 

didppm^ = dp{vPp)-pdpvl^, 



(3.36) 



where n is unite normal vector to dD, then 

/ j[-vgmk^w^v^ + 5idf,pv^]d^d^= J J [anm{w,p)y^]dsd^ 
l=-iD {=-iaD 

+ / IWgmkdxw''dxv"' + vdxgmkdxw^v"' - f>(9/3i;^]dxd^, 
i=-iD 

where anm(tu,p) = {-vgmkdpw^ +pSmfi)np. 

Recall the bilinear form and trilinear form on V{^) 



aQ{w,v) = J{vgmkdxw''dxv'^)dxd^, 

b{w, U,V) = / Bm{w, W)v"'^dxd^ = J gmk[dx{w^Vi^) + TT^w'u^v'^dxdC ^^'^^^ 



n 

■I 

n Q 
If the following boundary conditions are satisfied 



w\r^=o, onr, = r(urfaur|5=±i, 

<^n{w,p)\ri = h, OnFi = Tin U Tout, 



(3.38) 



then the variational formulation for (3.24) is given by 

' Findtu e V{n), p e Ll{n), such that 
aoiw,v) + {C{w,oj),v) + b{w,w,v) - {p,d„v") + {^^,v^) - uiire)-^agmk^,v"') 
+{d^^miw. e), + z/(P4.(e))9^«;^ - u{qmjie)w^ ,v"') 
{f,v)+<h,v>:^<F,v>, yveV{n) 

(3.39) 
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where 



4 The equations For the average velocity along the Rotating 
Direction 

We define the average along the rotating direction for the function ip[x^,x'^,$) in the coordinate 
(a;,^) in the domain 0, = D x [-1, 1] e E? 

1 

= \j ^{x,m ■■= ^- Vv'(x,0 e L2(fi) (4.1) 

-1 

It is well known that the divergence of a vector w can be written as under the coordinate 

^. _|_ dw^ _ 1 d{rw") ^ dw^ 

dx'^ r r dx"' ' 

Prom this it yields 

-1 

Since boundary conditions, 

h I = hi^'k=i - ^'l«=-i) = o,vw e v(i^), 

-1 

/ ^ "^dC = ^af-(rw«) = ^ + f := dTv.(w), 
-1 

where 

div2(t«) = ^ + — = -^(ru;"). (4.2) 

Therefore we assert 

M(divw) = div2(tu). (4.3) 

and the incompressibility becomes 



div2(a;\a;2w) = 0, (4.4) 
Taking into account the boundary conditions, 

i«|9+ua_U7tU76 = 0, -^|^=±i = -div2tu|{=±i = 0. (4.5) 

we get 

M(a^$(u;,e)) = 0, (4.6) 

Let make notation [w\ = w\^=x — w\^=-\,w = Mw. Then average equations of Navier-Stokes 
equations are given by 

diY2W = 0, 

-uE^iw) - vP^^iQY^ - vq^,{Q)w^ + 5iVpp (4.7) 
+Cm{w,uj) + MiBmiw, w)) = M{f)m + u{rs)-^agma[^] - Slip], 
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Let wv = axcrW^v"' + w^v^, w = w — w, then, it is clear that 

M{w^ - w^) = Mw = 0, M{ww) = 0, 

Hence 



Miw^w") = w^w" + M{{w^)w''), 



thus we conclude 
Finally, by virtue of 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



M{ww) = M(w{w + w)) = M{ww), 
it yields the reduced Navier-Stokes equations 

div2uJ = 0. 
We define the Sobolev spaces 

v{VL) = {u e h}{n), u = 0, onTtVb u r+r_, }, 

V{D) = {u e u = 0, on7o, see(3.13)}, 

By a similar manner as (3.39) the variational formulation for the reduced Navier-Stokes equations 
(4.11) is given as 



'Find w&V{D), p e L'^{D) such that 
aoiw,v) + iCiw,u;),v) - v(P^j{Q)dpwi + q^jw\v^) + h{w,w,v) - (p,5„v") 
= {-g^kM{dx{ww'') + T^k^ijiw'wi)),^^) 
_ +{u{re)-^agrna[%-]-Ssm\p\,v"') + {Mfm,v"'), Vi;ey(£>) 

[ (div2«;,g) = 0, VqG L\D), 



(4.12) 



where 



an{u,v) = (i'g.mkd\u'',dxv'^) = J vgmkdxu^dxv''^d-K, 
b{u,w,v) = {g^k{dx{u^w'^)+nheyw^),v'^), 



(4.13) 



5 The Equations for the Gateaux Derivative of the solu- 
tions of NSE with Respect to the Shape of Boundary 

In this section we consider the derivatives of the solution of NSE with respective to two dimensional 
manifold which is a portion of the solid boundary of the flow in the channel in turbo-machinery. 

Theorem 5.1. Assume that Surface 3 is smooth enough, for example, O G C^{D), then there 
exists a Gateaux derivatives (w := := of the solutions {w,p) of Navier-Stokes equations 

(3.24-) with respect to satisfy the following linearized Navier-Stokes equations : 



dimn := |^ 



r ' 



-v^w^ - „{re)-^a^ - vPf{Q)^ - 2^s-'e^gg, - vPf{&)^ - ,yq^{e)& (5.1) 
+g^fdpp + g^^d^p + C''{w,Lo) + iV«(«;, w) -|- N'^{w, w) -\- R''{w,p, 6) = 0, 
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w = o, on r^n{c = $fe}, 

'^w^ -pn = 0, on Tin n Tout, 



where 



R''{w,p,Q)r] := -2iy{re) "^OaVc 



Dg 



DO ^^iP 



Proof: The proof sec Appendix. 

Associated variational formulation for (4.3) is given by 

Find weV{n), peL^{VL) such that Vw e y(l^) 

ao(w, v) + (C(w, uj), v) + (L(«), 6), v) + h{w, w, v) - {p, d^v'^) + (^^f^, v) 



dw 



where b{-,-,-) and T(-, •) are respectively defined by (4.17) and (3.43). 



(5.2) 



(5.3) 



□ 



(5.4) 



6 2D-3C Navier-Stoke Equations on the 2D manifold 

As mentioned previously, for any ^ = const , there will correspond to a two dimensional surface S^. 
On the other hand, the three components of coordinate (x, ^) represent different meaning, the first 
two components x" are variables on the tangent plane to the surface 3^, which describe the flow 
direction in the channel, and the third component ^ is transverse variable which describe transverse 
flow through different manifolds. Therefor the Navier-Stokes equations (3.24) can be decomposed 
into two parts, the first is the operator on the tangent plane to the surface 3^ , which will be named 
"Membrane Operator", meanwhile the second is the operator along the transverse direction, 
which is named "Bending Operator" . Proceeding from this thinking, under the new coordinate 
the Navier-Stokes equations (3.24) can be rewritten as 



1 d{rw° 



J\f'(w,p,e) := -vAv 



div2W + 



9? 



0, 



V/3P + C'{w,u;)-vP{w,e) 



where 



+f^{riw,p,G))+B^{w,w) = r, 



B'-{w, w) — dp{w^w") + TTijW^w^ , tp\'>J^,P, 0) = w^w^ + rfp — uP^{w, 6), 



rf = [re) ^a, 



(6.1) 



(6.2) 



Let restrict the Navier-Stokes equations (6.2) on the surface SS^j. and adopt the Euler center 
difference to instead of the derivative with respective to ^ appearing in the bending operator A/j. 
Then we introduce several abbreviation about jump operator and finite difference operators. 



w{k) := i«|{={fe, [w]k := w{k + 1) - w{k - 1), 

(or)[tu]fe := w{k + 1) - w{k), {ov)[w]k := w{k) - w{k - 1), 



[w]k, T = ^fc+i - ^k, 



(6.3) 



and the corresponding different quotient represent as 



= ^k. = ^(«'"l«=«.+i -2t«"|£=€, +«^"l€=«._J 



(6.4) 
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Under this notations, we get 

-\ar[w'] - ^e^^{k - 1) + ^qijwHk - 1) - ^rj'p{k - 1) - ^w^ik)w'ik - 1) 

= K5j - '^q^w^k) + ^e^f^(fc) + + iu;3(fc)u;*(fc) + - 1), 

Ri{k - 1) = -ia,[«;*] - ^,Q^^{k - 1) + - 1) - ^TyM^ - 1) - ^w^{k)w'{k - 1), 

where Ur = 



So we finally conclude that. 

Theorem 6.1. The 2D-3C Navier-Stokes problem restricted on a smooth 2D surface is given 
by 

J\fi{k) := -uAw'ik) + 4(fc) + C'{w{k), w) + (p'l^V^pik) + \r}'p{k) 

+B^{w{k),w{k)) + \w^{k)w\k) = Fi.{k), (6.5) 
div2{w{k)) = —dl{w'^), {where div^w := ^da(rw°') , dr{w) = df.{w^)), 

with boundary conditions 

Crn(«^,P)l7.„ = hin, lin = U = ^fc} (6.6) 

a-n{w,p)\j^^^ = hout, lout = Tout U = ^fc} 



(6.7) 



where _ 

f a = l + r2|Ve|2, an{w,p) = -{v^-pn'')e^-{v%)e:i, 
L\{k) = {ar^l - -^qljjw^ik) - vr{w{k),Q) + f^e^9^«/(fc) 

= {-^P'l + Ti^^^Mw-^ + i-'^q^iQ) + o^rS] - ^qyw^, 

[ F^{k) = f^{k) + R^{k-l), 
Remark 6.1. There exists a second order differential operator in d^tp'^{'w,p, 6), 

and the term arw'^{k) in (6.5) is obtained by using the different quotient of second order. 

By a similar manner with (3.36), the equation satisfied by the covariant components of the 
Navier-Stokes equations are given as 

Ai{k) := QimJ^'^ik) = -ugimAw"'{k) + gimL^ik) + gimC"'{w{k),uj) + g„n4>'^PV pp{k) 
+\gimr]"'p{k)+gimB^{w{k),w{k)) + ^9imw''{k)w^{k) = gimF^{k), 

Simple calculation shows 

9im4''^^^ pp{k)v' = gimiS"'^ - S^^s-^Sl^-e,)d0p{ky = {gip - gise-^Qp)d0p{ky 
= {gap - g3a£~'^Qp)di3pv" + {gap - g33£~^Q0)d0p{k)v^ = Sapdppv"' 

= da{pv'^) -pdaV'^, 

gimiirv' = {-e-'^@aacxi3 + {rs)''^ ag3p)v^ + {-e-'^Q^gaa + {rey^agsav^ = v^, 
gim^iV^v' = dx{gimd\w'^v') - gimdxW^dxv' - dxgimdxW^v', 



Ai{k)v' = -iydx{gimdxw"'v') + vg.mdxw'^dxv'' + {vdxgimdxw"" + gimL':p'{k))v' 

+giraC"'{w{k),uy + da,{pv'') - pdo^v" +pv^ (6.8) 
+gimB^{w{k),w{k)y + ^girr,w\k)w^{k). 
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By using Green's formula and boundary conditions, we get 

J[-iydx{g^jdxiw^)vn + d^{pv'^)]d^= J (-i/gy^vj + pn„v«)ds 

D 7inU7out 

where n is normal vector to 71 = 7j„ U 7out, i.e., 

J n = n"ea + 0e3, n^ = ns = 0, 
\ <Jn{w,p) = -v^ +pn. 

Hence the variational formulation associated with (6.5) and (6.6) are expressed as 

' ¥mdw{k) &V{D), pe L^{D), such that 
ao{w,v) + mk),v) + {C{k),v) + b{w,w,v) + {p{k),v^ - d^v") 

= iGr{k),v), yvev(D). 

^ {(iiY2W,q) = {driw),q), VgeL2(L»), 



(6.9) 



(6.10) 



(6.12) 



where 

ao{w,v) = J{{iyg,jVxw'^xv^)dx 

____ 
= / fiaa^V xW^V xv^^ + r'^sQisiVxw'^^Vxv^ + Vxw^Vxv'^) + r'^ e^V xw^'^ xv^]dx, 

{C{k),v) = J{gijC'{k)v^dx = j2ioj[{Q^w^ - 602^{w, e))v'^ + ew^v^Jdx, 

{L{w{k), e)^v) = j[{ydxgimdxw^ + gimL^{k)y\AK ^^^^^ 

b{w,u,v) = {B{w,u),v) = (g,jB\w{k),w{k)) + gijlw^{k)w\k),vi) 
= {{dxiw^W^) + ao^pyl^ + er^e^v^) 

+{{dx{tu^w^) + n^j^w;'"?^;'^), er2n(t>, 9)), 
t {Gr{k),v) = {f,{k),v)+<ar,{w,p),h>\^, (by(6.7)) 

where 71 = ^in^lout, 

^=g7r + g^, diV2W = - '3^^ , 
V{D):={v\v€ h}{D), u = 0on7j, 

and 

Lij(e) = Urgij - ^gimq^j - l^girnQj^, (6.13) 

n?,(e) = 7rg.(e) + i%<5„„ nf^(e) = 4(6) + ^SsiSsj 
7 Pressure Correction Equation on the Blade Surface 

Noting that we must give value of [p] in the source term FT-(k) of equations (6.12), so the pressure 
on the surface must be supplied. Therefore, we recall the Navier-Stokes equations in invariant 

form 

' -lyg^'^VjVkW^ + w^VjW^ + 2s^^''gjmgkioj'"w^ + g'^Vjp = f\ , 

Let take divergence V, for (6.1)i and apply the identity Vfegjj = Vfcg*-' = Vfec'-'™ = 0, then we 
have 

-z/gj'^ViVjVfeu;' + div((w • V)tu)) + div(2w x t«) + ff'^ViV^p = div/. (7.2) 
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Because the Riemann curvature tensor vanishes in Euchdean space R^, therefore by exchanging 
the order of covariant derivatives, we have 



In addition, a simple calculation shows that 

div((w; • V)w) = Vi{w^SIjW^) = ViW^WjW^ + w^SJjWiW^ = ViW^VjW^ 



(7.3) 



(7.4) 



Prom (2,7) we have. 



C = 2u;xw = C'e„ = 0, = -2wrn(u;, 9), = 2u}e-\rQ2U{w,e) + '^), 
= -2ijr-H2rIl{w,e)+r^lUiw,e)) = -^-^l{r^IL{w,e)). 



On the other hand, the Laplacc-Bctrami operator can be expressed as 

Ap = ^d.ig^'^^dkP) = {re)-^[d^{rsg^f>d0p) + {d^{reg^^d^p) + dx{reg^^d^p)) 

+d^{reg'-^^d^p)], 

By using (2.5), we claim that 

= ^AS^ire^) - (2re.)gf^ - (6^ + rAQff^ + {re)-'a^], 
Assume that centrifugation force is the only exterior force, that is 

/ = —u} X u; X R = —LO^R, 

then 



(7.5) 



(7.6) 



On the other hand, we have 



R = rer = r{e2 - s "^8263), r 
div(/) 



r, r" = 0, r"* = —e ^r<d2 



1 _ n ».3 
re2n(w;,e) = l + + l = 2. 



-2 a; 



Summing up the above conclusions, we get 



US^ire^) - Mi^re,^ + (92 + rA9)p) + (r£)-ia|f ] 
-^|,(r2n(«,,9)) = -2|a,|2. 



(7.7) 



Next we consider the restriction of equation (7.7) on any surface 9^^. Noting that the Laplace- 
Betrami operator of pressure p in the new curvilinear coordinate system (.x" , ^) can be split as the 
sum of two operators , membrane operator on tangent space and the bending operator along the 
rotational direction 



-Ap = -Amp - Afep, 

-AmP = - — g^(r£g^) = --g^(rg^), ^ 
-A,p = -i[(re)-ia|f + |(2r8,^ + (92 + rA9)p)]. 



(7.8) 
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We approximate the derivatives with respect to rotational variable in (7.7) by the difference quo- 
tients defined by (6.5), and then restricted it on the 9^^., finally we get 



dnP = fn, other boundaries . 



(7.9) 



where 



Mr) = -2re,dl{^) (92 +rAe)4(p) + ^4(p) 
-2ujr-'-§-^{r'U{w{k),Q))-2\co\' 



Ysw^ := dliw^) - reS2xIl{w{k), 9), 
Vaw^ := dliw^) + r-^w'^{k) + r92n(w(fc), 9). 

Specially, we consider the restriction of (7.7) on the surface 3±i. When ^ = ±1, we have 

t('l5=±i = 0, daW = 0, VqW* = 0, 

V3U'3(±1) = ^(±1) = -{daw°' + r-iw;2)(±l) = 0, 

n(u;,9) = 5^i^=0,(see (3.19)). 



(7.10) 



At present, equation (7.9) becomes 



13 (r ^'P ) 

= (2^0A^ + (02 + rA9)p) + (re)-ia|f - 2{lo)\ 



(7.11) 



Furthermore, we replace the derivatives ^ by difference quotient and apply the boundary condition 
w|5=±i = 0, then 



^|{=-l = i(t«|£=-l+r - i«|«=-i) = ^(w|£=-l+r), 
^|i=-l = ^(t«|5=_i+2r - 2«,|j=_i+,), 



(7.12) 



Borrowing the notations in (6.3), when ^ — ±1, 

r d^,{w) = ^{-l)^lw{-l + T), (u;(-l) = 0), 



1 dl,{w) 
Similarly, when ^ = 



9C 
a{2 



^ "^-(-1) = U-w{-l + 2t) - 2ti;(-l +t)), {w{-l) = 0), 



\ dl{w) = ^(1) = U'^il - 2t) - 2ti,(l - r)), {w{l) = 0), 
In addition, Owing to ^ = — 1 and ^ = 1 are both sides of blade , hence 

p(-l-T)c.p(l), p(l + r) =p(-l). 
therefore, we rewrite (7.13), (7.14) as 

f d_i(p) := ||(-1) = i(p(-l + r)-K-l)), 

dlM ■■= |f (-1) = ^(P(-1 + r) - 2M-1) +P(-1 - r)) 
= ^{p{-l + T)-2p{-l)+p{l)), (p(-l-r)=p(l)), 

I rf2,(p):= J,(p(-l + r)+Kl)), (M-l-r)=p(l)), 



(7.13) 



(7.14) 



(7.15) 



(7.16) 
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^(l) = i(p(l)-p(l-r)), 



Mp) — w 

dlip) ■■= |f (1) - Mpi^ + ^) - 2P(1) - r)) 



(7.17) 



Summing up and introducing p_ = p{—l), from (7.5) the equation of the pressure on the ^ = — 1 
surface is give by 



1 9 f^dp- 



(r^) + arP- = /-(t), 

/_(t) = -2reAd-i(|S,) - (62 +rAe)d_i(p) + f^dl,{p) 
+d^i(u;3)-2|w|2, 

By a similar manner we can obtain the equation of pressure on the surface ^ = 1 

-igf^(r§|i) + = /+(t), 

U{t) - -2reAdi(^) - (62 + rAe)di(p) + ^d?(p) 
+d?(u;3)-2|cj|2 



(7.18) 



(7.19) 



Next, in order to inspect the reliability of the method, we make some numerical simulations 
of pressure field by using the pressure correction equations on the blade surface (7.18), (7.19), and 
boundary conditions in (7.9). The low speed large-scale centrifugal impeller of The NASA is used 
as the example ([22]), and some comparison with FLUENT's conclusions is diagramed as below, 



P:pa 



FLUENT6.3 



P:pa 




0.1 0.2 0.3 0.4 




0.1 0.2 0.3 0.4 



(a)Pressure distribution on the pressure surface 
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(b) Pressure distribution on the suction surface 
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P:pa 



FLUENT6.3 P:pa 



LB-2D 




(c)Pressure distribution on the median surface 
Fig5. Numerical comparison of the Fhient's conclusions and the Pressure Correction method 
[Numerical results was completed by Chen Hao, Energy and Power Engineering College of Xi'an 

Jiaotong University] 

Where the LB-2D on the RHS are the results from the Pressure Correction method. 
Fluid power, Fluent evaluates to 13973 watt, our method is 13975 watt. 

8 Steam Layer in Domain decomposition and the Bi-parallel 
Algorithm 

Next we we consider the decomposition of the flow passage, 

o = 7? X {-1,1} = ^{i? X = J^f^fc 

k k 

where — l = ^o<^i<'''< (.tn = 1- In the subsequent ftk is called "stream layer". 




5o 



Fig 3: decomposition of the flow passage and angular expansion 

The new method is to solve the 2D problem (6.10) with respect to the velocity and the pressure 
(w^p) and the pressure correction equation (7.9) on the interface 3^ of two stream layer flk and 
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Qfe + 1, which is a 2D manifold. The general parallel algorithm can be used to solve these 2D- 
problems (6.10) and (7.9). 

Bi-Parallel Algorithm: The Bi-Parallel Algorithm means that we adopt the parallel algorithm 
to solve the problems (6.10) and (7.9) along the two directions, i.e., on the 2D manifold 9^ and 
along the direction ^. On the specified 2D manifold 9^, the general domain decomposition method 
or data parallel algorithms can be used to implement the parallel algorithm. On the other hand, the 
problems on the 2D manifold corresponding to different discrete parameter ^i.,k = 1, • • • , m, 
can be solved at the same time, which forms another parallel. 

As new method is applied to solve 3D-viscous flow in turbo-machinery, all interfaces Sj^^. have 
the same geometry properties, i.e., the same ^a/3) * • ' . On the other hand , when this methods is 
applied to other 3D-flow, for example, circulation flow through the aircraft, geophysical flow around 
the earth, in this case the interface surface 9^^. have different geometry properties. Suppose the 
next interface Sj^.,.! is generated by a displacement t] of the previous interface , then the new 
fundamental forms {aafi{'n),bap{'n) can be computed by the following formulas. 

Theorem 8.1. Assume that Q is a smooth surface in 5R^, a^^, 6q,^ are metric and curvature tensors 

respectively. Given a smooth displacement field r] = rj^Ba + rpn of^, we get the new surface '^{rj), 
and use symbols aai3{r]),bai3{r]) to denote the metric tensor and the curvature tensors of the surface 
$S(r/). Then a simple calculation shows that they can be expressed as(see [3,4]), 



where 



bafsij]) 



Oafj + 2 Eafj (??), 

baP + PapiTl) + Q'ilii'n), 

{boc0 + Papir])){q{v) - 1) + qivMMvJdiv) + <^(??)TOa(r/) 

-iPaffM+KffVx Vn Va V% 



Eal3 (r?) = la^iv) + ^]aXa Va V V 13 Va V V/3 V ] 

* ,0 

Papiv) =VaV/3 V + baa V/3 , 

V^ € =V/3 - b^^v^ V^ r -- 



lapin) = |(a^A Va + «aA V/3 



(8.1) 



(8.2) 



Pa^iv) =VaV0 V'' - b^ V/3 ??^ 

* 

(t>a0{v) = bap + Pa0{v)+ ^aB^ \ V 



daiv) 
d{ri) 



do{7j) = l + d{Tj), 



7o(?7) + det{\' 



mi{r]) = 2 





Vi rf 





V2 ri 





m2{'n) = -2 



, 
Vi T 





, 
V2 V 





Vi n V2 T 

where (5^^(t?) is a remainder term which order is higher than 1. 
Proof: The proof is omitted. 



Vi 11 V2 n 



1.3) 



(8.4) 



□ 



9 Existence of Solution of the 2D-3C Variational Problem 

In this section we study the 2D-3C variational problem (6.10) on the manifold First, let 
dD = 7^ u 7o and introduce the Sobolev space V{D) deflned by 

V{D) = {w\,w = {w", w^) e H^{D) X H\D) x H\D), w\^^ = 0}, (9.1) 
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equipped with the usual Sobolev norms 

r HiD = Ej:\\daw^\ij„ ii«'ii^,z. = Eikn^,z, = E/ki'dx, 

I ct 3 i i D (9.2) 

I \\hId = \mId + m\Id- 

It is clear that the variational problem (6.10) is a saddle point problem. In order to regularize 
it we introduce the artificial viscosity r] such that 



' Findw e V{D), p e L^{D), such that _ 
{arW, v) + ao{w, v) + {L{w, Q),v) + {C{w, oj),v) + b{w, w, v) — {p, divv) 
=< Gr,v >, \fv e V{D), 
[ V{P, q) + (div2W^ - dr, q) =0, Wqe L^D), 



(9.3) 



where 

< Gr, V >= if^, v)+ <h,v>\^„ (9.4) 
Obviously, problem (9.3) is equivalent to 

Findw e V{D), such that 

Ao{w,v) +r]~''^ {div2W, divv)) + b{w,w,v) =< G,v >, \/v € V{D), (9.5) 
p = r]~^ [-div2 w + dr{w)], 



(9.6) 



where 

Ao{w, v) = {arW, v) + ao{w, v) + {L{w, 6), v) + {C{w, w), v), 
<G,v> = {f^,v)+ < hin,v > |r,„ + (r?"MT-,divt;)- < r]-'^dr,v >, 

Our first objective is to show that the bilinear form Ao{-, ■) defined by (9.6) is y(I?)-elliptic. 
Theorem 9.1. Let D be a bounded domain in B?, the injective mapping ^{x) defined by (2.1) 
satisfies 3? e C^(D), |-D'3fJ|oo,D < ko, and the two vectors = dorSi- are linearly independent at 

all points of D- Let 70 be a d') -measurable subset of ^ = dD and d'jo > 0. Then there exist the 
constants C(fco), C'o(fco) depend upon x such that the following equality and inequalities hold 

(i) aoiw,v) — aa{v,w), Ww,ve V{D), 

(ii) \ao{w,v)\<iyil + rfkl)\w\i,D\v\i,D, 'iw,veViD), (9.7) 
(Hi) \ao{w,w)\ > lylwll ^, yw,v& V{D), 

( (iv) \Aoiw,v)\< C{ko)\\w\\i^D\\v\\o,D, yw,veV{D), 

(v) Aoiw,w)> {:^:^-Coiko)ko)\\w\\ljj + {iy-Ci{ko)ko)\w\ljj ^g_g^ 
+ / ^ar(w««;")dx, \/w e V(D), 

D 

Proof: Firstly, from (6.11), we have 

ao{w,v) = J[iygijVxw'Vxv^]dx = /[z/a^^VAW^VAV^ 

+{erVxw^)irei3yxv'^) + {erV^xv^){rQi3Vxw'^) + {er'V xw^){re'V xv^)]dK 
= yJiScDVxw'^Vxvl' + irQc,Vxw'')ir&^Vxv^) 

d _ _ „ „ „ „ 

+{rQ0Vxw'^){reWxv^) + {reV xw^){rQpy xv^) + {reW xw'^){re{V xv^)]d^, 
By combining we have 

ao(iu, v) = v j [5apVxw"Vxv'^ + (re«VAw" + reV xw^){rQpV xv^ + reVxv^)]dl^ (9.9) 
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hence, we get 

aoiw,v) = ao{v,w), \ao{w,v)\ < + rlkl)\w\i,D\v\i,D- 



Then (9.7) is proved. Next wo consider (9.8). Obviously {iv) is vahd, thus we just need to prove 
(v). Indeed from (3.26), (3.27) and (3.35) we assert 

f L^{w, 9) = {r-'S2/3Sc.a + r^QaOafi - r^ea^fia)df!w'' - (2ree^a2a + er^ef!„)d0W^ 
L3{w, 6) = r'^eQc^fsdpw"' - 2e'^r^e2Q,3daw^ + (er^QaqZ + r"^ s'^ q^)w'^ , 
aaaQp + r^eqj = r-^62i3S2a + S2isQ2Qa + 2a2<x(<52/3|Ve|2 - aGae^s 
aaaQa + r'^eql = ra2„Ae + £aQ2{rQ2 - 2a2a), 

{er^Qcq^ + r^elql) = s{e2S20 + 6^3) + 2^82/3 + r2£e2(2(52;3|Ve|2 - 0626^) 

+£r3e2e;3Ae - r^ee2ex&xp, 

{ (er^e^gf + r'^e'^ql) = r^e^{rAe - 063), 
therefore, 

Lj{w, e)w^ = [{r-^52fiSaa + r^QaOafi " r^Qa^ i}a)djw'^ - (2^6^02^ + er^ep„)dpw^]w'' 
+[{r-^S20 62a + ^2^626^ + 2a2a(<52/3|Ve|2 - aOse^s 
+re^Ae - rexQxfi) + 00220,6^ + 2re,e20 + r^eadpAe)w^ 
+{ra2aAe + £ae2(re<^ - 2a2a))w^]w'^ 
+[r^eea)3d,3W°' - 2e^r^e2@f!d)3W^]w^ 
+[{£{02620 + Qp) + 2reQ20 + r2ee2(2(52/3|Ve|2 - 0626^) 
+er3e2e^Ae - r^sQ2Q\O\0)w^ + {r'^e^{rA@ - a@2))w^]w^ 

By simpHfying we get 

Lj{w,Q)w^ = r-^^{w°'w") +r-'^w'^w'^ + P^{w,Q)d0W°' + Pi{w,e)d0W^ 



+(5a/3(e)w;^«;" + g3/3(e)u;'^W;3 + Q33(e)u;3^3 



where 



r Pi{w,e) =r^e&^fsw^ + r^Q„e^p-e„ef,,)uf 



pp{w, e = -£[2r3ee2e^w3 +jy.2(e^^ + ree2cre/3) 

Qa/3(e) = d20{Q2Qa + 2a2a|Ve|2) + ae;3(a22e« - 2a2„e2j 

+ (2r£52Ae„ + 2^12.6^)6^^ + re^(2a2„ + re„)Ae, (9.11) 
Q3/3(e) = s[r{2a20 - 620) AO - r^e2exexfi + Srae2&0 + Op 

-((1 +Jl)(52/3 + 002/3)62], 

I (333(0) = £V(rAe - 062), 

From the assmuptions of the Lemma, we clare that there exist two constants Cj(fco),i = 0,1 
independent of it;, such that 

{L{w,Q),w) > y i5,(w«w;")dx-Co(ko)ko||i«||o,D-Ci(ko)ko|i«|?,D 



In addition 



{C{w, u)),w) = 2 J {^ ^ w)wdx = 0, 

D 
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Then, from (9.6) and Ur = tj^^ > ^^"2^.2 , we get 

Ao{w,w) — {arW,w) + ao{w,w) + {L{w, 6), to) + {C{w,u}),w) 

> (^w - Co{ko)ko)\\w\\l o + Ci{ko)ko)\w\l ^ + / ia,(«;«to«)dx, 

^ D 

the proof is ended. □ 

Lemma 9.1. Under the assum.ptions in Lemma 9.1, the trilinear form bo{-, ■, •) is continuous, i.e., 
there exists a constant M{0,D) independent ofw,u,v, such that 

\biw,u,v)\ < M\\w\\^.^^^\\v\\^.^^^\\n\U,D,yw,u,v € V{D) (9.12) 



Proof: Thanks to the Holder inequahty 

\w^Vxu''vf\^d^ < ||t«||L*(D)||v^||L^(D)||Vu"||o,D 

and the Sobolev embedding theorems 

||«||l4(d) < C\\u\\^.^^^, Mlh^,) < CM^.^^y (9.13) 

Further, by using tlic Cauchy's inequahty, we derived the conclusion. The proof is completed. □ 
Remark 9.1. It is clear that we have 

aafiW^ Va w^w^ =Va {aa0W^U°'V^) - aafjw'^ Va (w^w^) 
=div {\w\w) — \w\ div w — actpw"w^ Va w", 

Hence 

* 1 * 1 * 

o-apw^ Va w"w^ = - div {\w\w) — divw. 
Considering (6.6) and the boundary conditions of element in V{D), we claim 

Mwa,wo,wo)\ < C(||«;o||i4(B)|| divwWlj, + ||t«o|li3(^^)), (9.14) 
where the Gauss theorem is used. 

Theorem 9.2. Under the assumptions in Lemma 9.1, for the given (G, rfg) S V*{D) x H~^(D), 
if F satisfies the following condition, 

\\F\\*<'mus, with < F,v >=<G,v > -r]-'^{dl,divv), (9-15) 
then there exists one solution of the variational problem (6.10) which satisfies 



Furthermore, if 

m* < (9-17) 

then problem (6.10) has a unique solution in V{D). 



23 



Proof: We begin with constructing a sequence of approximate solutions by Galcrkin's method. 
Since the space V{D) is separable, there exists a sequence {ip^,m > 1) in V{D) such that: 1). 
for all m > 1, the elements ¥>i, • • • are linearly independent; 2). the finite linear combinations 
J^^ifi dense in V{D). Such a sequence {(p^,m > 1) is called a basis of the separable space 

V{D). 

Next we use Vm to denote the subspace of V{D) spanned by finite sequence cpi, ■ ■ ■ , tp^. Then, 
we can construce the approximating problem, 

f Find Wm e Vm such that , , 

\ Ao{Wm,v) + bo{Wm,Wm,v) =< F,V >, \/ V € Vm- 

If we set 

m 

Wm = ^QVSj, 
i=l 

then wc find that problem (9.18) amounts to solve a system of m nonlinear equations with m 
unknowns q. For each m problem (9.18) has at least one solution. Indeed, we can introduce the 
mapping Mm ■ Vm ^ Vm, 

(A1,„(m),(PJ = Aoiu,ip^) + bo{u,u,cp^)- < F,ip^>, l<i<m, 

where (•, •) is the scalar product in V. Hence, Wm G l^n is a solution of problem (9.18) if only if 
Mm{wm) = 0. Since 

{Mm{u),u) = AQ{u,u)+bo{u,u,u)- < F,u>, VuG Vm, 

it follows that 

iMm{u),u) > (^||w|li,i, - M\\u\\l^ - ||F|U)||u||i,^ (9.19) 

Hence, if choosing 

^ ~ MC V ^MC' M ' 
then we get for all u € Vm with ||tx||i,£) = p, 

{Mm{u),u)>Q. 

Moreover, M.m is continuous in Vm, and the space Vm is finite dimensional, we can apply Corollary 

1.1 in [17], there exists at least one solution Wm € Vm of problem (9.18). 
Furthermore, we have for any solution Wm to (9.18) 

= {Mm{Wm),Wm) > (^ ||l«m||l,£> - ||l«m " ll-P'll*) ||l«m|| 

therefore, 

^\\wmKD-M\\wm\\ln-\\F\U<0- 
It follows that, when denoting by y = Htfrnlli,!)) 



MC' -^MC' M ^ MC - y^MC' M ' 
i.e., 
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This shows that the sequence {w„i) is uniformly bounded in V. Therefore, we can extract a 
subsequence, still denoted by Wm), such that 

Wm (weak) in V{D) as m ^ +00. 

Then, the compactness of the embedding of V{D) into L'^{D)^ implies that 

■J^m (strong) in asm ^ +00, 

the remainder is to prove b{-, ■, •) is weakly sequence continuous, i.e., bo{wm,, Wm, v) — )■ 6o(ii'», tu*, f). 
To do this, we recall 

V — {u £ C°^{D) satisfy the boundary conditon (3.36) } 
is dense in V{D) and 

bo{wm, Wm, v)= I aaffW^v^ ax^w^rf dl - bo{Wm, V, Wm)- 

For any v G V, then v € L°°(D) L°°(7i), dx^'V^ € L°°{D), and the convergence relations 
lim WmWm = w^wl are satisfied in L^{D) and ^^(71) respectively, therefore 



m— >oo 



lim bo{wm,Wm,v)= L aa8W^v<^axaW^n'^dl - bo{w*,v,w*) = bo{w^,Wt,v), \/v € V. 

m— >cx) '1 

Next for all v G V{D), by virtue of the density of V, and taking the limitation of both sides of 
(9.18) implies 

Ao{w^,v) + ba{w^,w^,v) =< F,v >, Vt; e y(£>), (9.21) 

that means is a solution of problem (6.10). 

In order to prove (9.16), by a similar manner we get 



MC' - V ^MC^ M ' 



l|w*||i,D> T77^ + \/(T77^)2--r^,or < - , / ( )2 



MC y-MC' M ' " ""'"^ - MC y^MC" M ' 

Obviously, the first one of the above inequalities is contract to (9.20), thus only the second one is 
true, that is (9.16). 

Next we prove the uniqueness. In fact, if there exist two solutions and of (6.10). Let 
e* — — Wt,, then 

Ao{e^, e*) + 6o(e,, w*,e^) + bo{w^, e,, e*) = 0. 
Owing to condition satisfied by tw* and u;, and (6.17), we get 

>(^ - 2M(^ - JA^ - &)lle f -2J(^)2-M£ 

This yields ||e*||i,£) = 0. Therefor, the solutions is unique. The proof is complete. □ 

Theorem 9.3. Let {wo,po) and {Wjj,prj) be the solutions of (6.10) and (9.3), respectively. If F 
and = sup \H\ satisfy the condition 

D 

Co - 2Mp - 2r)-^Hl > C2 > 0, (9.22) 
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then the following estimates are valid 

\\wo - Wr,\\i,D + \\po - Pri\\o,D < max(C3, €4)11, (9.23) 

where 



C + 2Mp.. ,, ^ (C + 2Mp)2 

llPollo.D, W - 

and I3{) is the constant in the inf — sup condition. 



Cs='^^^\\po\\o,n, C,= '-Z.:2 \\Po\\o,n. (9.24) 



Proof: From the assumption wo have 

' Findioo € V{D),po e i^(-D), sueh that 

< ao{wo,v) - {po,divv) + bo{wo,wo,v) + {l{wo),v) =< G,v >, € V{D), (9.25) 
, (div Wo - 2Hwl + dl,q) = 0, Vg € L^(-D), 

and 

' Findu?,, e V{D),p,-, e L'^{D), such that, 

< ao{wn,v) ~ {prp div v) + bQ{wn,Wn,v) + {l{wn),v) =< G,v >, G V{D) (9.26) 
. T]{pr„q) + {divWr,-2Hw^ + dlq)=0.yqe L^D). 

Next, we denote e* = — = po— Pn- subtracting (9.25) from (9.26)then yields 

* 

ao{e^,v) + bo{e^,wo,v) + bo{wr„e*,v)- {s^,divv) =0,Wv G V{D), (9 27) 

(div e,,q) + {-2Hel, q) + v{s.,q) - v(Po,q) = 0,^ q ^ L'^{D). 

choosing = e*, = in (9.27) and summing the above two equations, we obtain 

oo(e*,e*) + 6o(e*,tuo,e*) + bo{wj,,e^,e^) +r]{s»,s^) - r]{Po,s») - {2Hel,s^) = 0. (9.28) 

Noting that (9.7), (9.8), and (9.12), we have 

(Co - 2Mp)||e*||^_^ +r?||s*||^ j3 < (r?|bo||o,r> + 2sup |F|||e*||o,D)||s*||o,D. (9.29) 

Furthermore by using Young's inequahty we get 

2sup l-ff |||e*||o,D||s*||o,D < o^lk*|lo,z? + 2»?-'H^||e*||?n, (9.30) 

D ^ 

therefore from (9.29), (9.30) 

(Co - 2Mp - 2r,-i7J2)||e,||2_^ + ir,|ls,||2_^ < r,||po||o,D||s*||o,D, 
Noting that the condition (9.22) are satisfied, hence 

II ||2 / |bo||o,D|| II /noi\ 

\\e*\\i D<V Jl \\s4o,D- (9.31) 
On the other hand, the inf — sup condition means that 

P0||S*||0,D < sup ' ||.^|| " 

vev(D) " 

< sup {\\v\\^]j\[ao{e^,v) + bo{e^,WQ,v) + bQ{wr,,e^,v)]\ (9.32) 
vev(D) 

< (C + 2Mp)||e*||i,D. 
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Finally from (9.31), (9.32) we have 

||e*||i,r( < C3?7^, ||s*||o,r> < (^4??, 

where 



C + 2Mp {C + 2Mpf 
^3 - ^ o — llPollo.D, w llPo||o,r>- 

Thus the theorem is proved □ 



10 Finite Element Approximation Based on Approximate 
Inertial Manifold 

In this section, we focus on the variational problem for the 2D-3C problem (6.10), called 2D-3C 
variational problem, 

f Find wq e y{U), such that , 

where 

Aq{wq,v) = ao{wo,v) + r/"i(div t(;o,div v) - r]~^{2Hw^,div v) + {loiwo),v) 

* * * 

= ao{wo, v) + 77~^(div wq, div v) + {u — r]~^){2HiUQ, div v) 

ao{wo, v) = 2iy{a°'^a>^''-fx^{wo),jc.p{v)) + y{a°'^ Va Wq, V/3 v^) 

bo{wo, Wo, v) = {aafjw^ Va - 26„/3Wo Wq, vI^) + {wl^ V/3 + ba/3W^w[l,v'-^), 
{l{wo), v) = (a«^P(u;o), v^) + {l^{wo),v^) 

* 

<Gr,,v> =<Fh,v>+ J^^ [anav" + (T„3V^]d-y + {aaffd^ml, v^) - r]-'^{dl, div v) 

(10.2) 

We now consider the finite element approximation of the 2D-3C variational problem (10.1). 
Assume that Vh and M/j are finite element subspaces of V{D) and L'^{D) respectively. Introduce 
the product space Yh = VhX M^, obviously which is a subspace of F = V{D) x L'^{D). 

Then the standard Galerkin finite element approximation of (10.1) is defined by 

r Yindwh e Vh, such that , . 

\ Ao{wh,v) + hQ{wh,Wh,v) = {Gh,v) "iv^Vh 

As usual, we make the following standard assumptions on the finite element subspace 

(HI) Approximation property 

inf {h\\u - Vh\\i,D + \\u - Vhh,D + h\\p - g/.||o,D} < C h^^^ {\\u\\k+i,D + \\p\\kd} 

for any {u,p) €Yn {H''+'^{nY x /f'=(f2)), l<k<l. 
(H2) Interpolation property 

\\v - ihvWi^D + \\q - Jhq\\o,D < C h''{\\v\\k+i^D + \\q\\k,D) 

for any {v, g) G F n (if'^+^(f2)'' x H''{Q)), I < k < I, where Ih, and Jh are some interpolation 
operators from V{D) and L'^{D) into Xh and M^, respectively. 
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(H3) Inverse inequality 

||t^^.||i,D<C/i-i||t;^||o,D ^Vh^Vn. 
(H4) {Vh,Mh) satisfies LBB-condition 

(q, divf ) 

mf sup T-- 7-- > ^ > 0, 

I^MnvaXn \\<lh,D \\v\\l,D 

where /3 is a constant independent of h. 

The fohowing optimal error estimates of the Galerkin finite element approximation are well- 
known (cf. [17]), 

Theorem 10.1. Suppose Wq G V{D) C\H^+^{Df is a nonsingular solution of (10.1) and the 
finite element subspace Vh satisfies assumptions (HI) ~ (HA). Then there exists a solution Wh 
satisfying (10.3) such that 

h\\wo - Wh\\l,D + - Wh\\o,D < Ch'^+'^\\wo\\k+i,D- (10.4) 

Next we try to improve the error estimation. To do that we rewrite (10.1) in operator form. 
Let T : V{D) — > V*{D) denote the 2D-3C Navier-Stokes operator on the manifold 9 via 

<T{wo),v>:= Ao{wo,v) + hQ{wo,wo,v)- <Gr,,v>, \/v€ V{D). 

It is obvious that T{wo) = is equivalent to (10.1). The operator form of finite element approxi- 
mation (10.3) is 

<Th{wh),v>:=Ao{wh,v) + bo{wh,Wh,v)-<Gn,v>, Vt; e Vh{D). 

Therefore, Fh{wh) = is equivalent to (10.3). Furthermore, it is easy to show that F{w(^) and 
^h{wh) are Frechet differentiable and the Frechet derivatives at Wo and Wh along direction u are 
given by, respectively 

Awoi'^^'") ■= {DT{wo)u,v) = Ao{u,v) + ho{u,WQ,v) + ho{wo,u,v), Vm, v e V{D), 
■Awh{u,v) := {DJ^h{wh)u,v) = Ao{u,v) +bo{u,Wh,v) + bo{wh,u,v), Vw,i; e Vh{D). 

It is well known that Wq is a nonsingular solution of (10.1) if and only if DJ^(wq) is an isomorphism 
on V{D), furthermore, equivalent to Awoi'i') satisfies the inf — sup condition(weak coerciveness) , 
i.e., 

mf sup ° ||,,|| — > ao > 0, mf sup n,,,, ° n,,,, — > qq > 0. fin 5) 

In this case, for any / e V*{D), the variational problem 

Find we ^(Z)) such that ^'1nf^^ 
Awo{u,v)=<f,v>, \fveV{D), ^^"-^^ 

has a one and only one solution. Similarly, Wh is a nonsingular solution of (10.3) if and only 
if DFh{wh) is an isomorphism on Vh{D), equivalent to Awhi'i') satisfies the inf — sup condi- 
tion(weak coerciveness) 

. Aw (U,V) . An, iU,V) 

uetin, viZm m.MviU^o > > , ^ -P^^ mu^EWKB ^ "'^ > (10.7) 
In this case, the variational problem 

r Findit^iS 14(D) such that 

\ AwA'^h,v)=<f,v>, VvhGVhiD), ^'^-^^ 
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has a unique solution Wh for any / e V^{D). Condition (10.5) is equivalent to 

\\DJ^{wo)\\civ,v)<ao'- (10.9) 

The next theorem shows the uniqueness condition to insure the finite element solution Wh of 
(10.3). 

Theorem 10.2. Assume that the assum,ptions (HI) ^ {HA) are valid, and wq is a nonsingular 
solution of (10.1). If the finite element mesh h is small enough such that 

2MCao^\\wo\\2,Dh < 1. (10.10) 

Then, solution Wh of the finite element approximation problem (10.3) is nonsingular. 

Proof: In fact, from the above explanation in this section, it is enough to prove that 

\\DTh{wh)\\ciVu,v^)<Po^- (10.11) 

Therefore, noting that 

£ : = \\DF{wo) - DFh(wh)\\c(v,v) 

_ '' bo{Wo-Wh,U,V)+bo{U,Wo-Wh,V) (10.12) 

< 2M||wo - Whh,D < 2MC\\wQ\\2,Dh. 
Set B = {DJ^{wo)}~'^{DJ^{wo) - DFh{wh)}. Then from (10.9) and (10.12), we have 
DTh{wh) = DJ^{wo){I - B), 

\\B\\c,y,y) < ao-^2Mq|«,o||2,../., ||(/-S)-^IU(v,^) < ,_,MC.-l^o||... ' 

\\DTh{Wh)\\c{Vy) < Wl-2MCao-\|,<,o||2,z,/^- 

substituting (10.10) into the above inequality, then we get DTh{wh) is an isomorphism on Vh, 
hence Wh is a nonsingular solution of (10.3). □ 

Theorem 10.2 shows if mesh size h is small enough, then we have 

mf sup — n— > -ao > 0, mf sup t--^ r-- > -ao > 0. 

(10.13) 

Next, assume Wh is a nonsingular solution (10.3). We define a projection Ph : V{D) — )• Vh{D),\lw G 
V{D) through 

Awdw - PhW,v) = 0,yv e Vh{D). (10.14) 

Since Wh is a nonsingular solution, then there exists a unique solution of (10.14). Consequently, 
V can be decomposed into the direct sum of two subspaces: 

V{D) = Vh{D)®Vh{D). 

This meas that for any w G V{D), we have 

W = PhW + Pj^W =Wp+ Wq, Wp e Vh{D), Wq € Vh{D). 

It is straightforward to show that 

( AwH{Wq,Vp) =0, \/VpeVh{D), Aw,,iw,Vp) = Aw,,{Wp,Vp), noT^^ 

I \\wqh,D < Ch''\\w\\k+i,D,yw e VnH'^+\D)\ ^'^-'^^ 
Next, we present some technical lemmas. 
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Lemma 10.1. there exists a constant independent ofu,v,w such that 

11 11 
|6o(w,w,'u)| < C||ti||^^^||w||i^_^||w||i,£)||t;||^_^||u||i^_^, \/u,w,ve V{D). (10.16) 

Proof:By virtue of Holder inequality we get 

\bo{u,w,v)\ < C||w||o,4,r>||w||i,D||i;||o,4,D, 

furthermore the Ladyzhenskaya inequality shows that 

MoA,D<C\\u\\l,^^\\u\\l^^^. 
Therefore we prove the lemma immediately. □ 
Since wq € V{D), it can be decomposed into wq = wop+woq with wop e Vh{D), w^q G Vh{D). 
Lemma 10.2. The following estimation is valid, 

2M 

\\Wop - Wh\\l,D < \\Wo - WhWl'jjWwo - Wh\\o°D, (10.17) 

ao 

where 

1, for the Homogenous Dirichelet B.C. = {B.C. I), 



£i = i 3 



2 ' 



for Mixed B.C. = (B.C. 1 1), 



_ J 1, for the Homogenous Dirichelet B.C. = {B.C.I), 
^° ~ I i ^«2;ed B.C. = {B.C. 1 1). 

Proof: Firstly, equation ((10.1) can be rewritten as 

■Awhii^Og,v)+Awu{'Wop,v) + bo{wQ-Wh,wo-Wh,v)-bo{wh,Wh,v) =< Gn,v>,yv e V{D), 

(10.18) 

meanwhile (10.3) as 

Aw,Awh,v) - bo{wh,Wh,v) ^< Gr„v >y V e Vh{D). (10.19) 
Let V G Vft, subtracting (10.18) from (10.19) and using (10.15) with w = Wq, we derive 

Awh{wop - Wh,v) = -bo{wo - Wh,wo - Wh,v). (10.20) 
Since Wh is nonsingular, (10.13) shows 

irv^ < 1 Aw^iWop-WH,v) _ 1 „,,^ -boiWo-WH,Wo-WH,v) 

that is, 

\ao\\wp-w,\U,,< sup l^oiwo-w^w w,,v)l _ 
veVh{D) " 

By using (10.16), for any v G Vh{D) we get 

\bo{wo - Wh,wo - Wh,v)\ = \bo{w„ - wi,, v,wo - Wh)\ 

< M\\wo - Wh\\o,D\\wo - '«'/i||i,D||'y||i,D, for B.C./, 
\bo{wn - Wh,wa - Wh,v)\ = \bo{wo - Wh,v,Wo - Wh)\ 

< M\\wo - t«;,,|||_^||u;o - Wh\\lD\\v\\i,D, iovB.C.II, 
Summing up the above relations we draw the conclusion (10.17). □ 



(10.21) 



Next, let the mapping 0(-) : Vh{D) — )■ Vh{D). we define the manifold M as the the graph of a 
function (p, that is, A4 = Graph^, then problem (10.1) can be rewriten as 

I Fmd(p{w) e Vh{D), such that ^ 22) 

We first give an approximation property of the solution of (10.22). 
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Theorem 10.3. Suppose that the finite element spaee Vh satisfies the assumptions (HI) ~ (^4). 
Then there exists a mapping <p{w) defined by (10.22) which is a Lipschitz continuous function with 
the Lipschitz constant I = l{p), and </>(•) attracts any solution wp, of (10.1), i.e., 

{H5) \\(t){Wi)-ct){W2)\\l,D<l{p)\\wi-W2\\l.D. 'iwi,W2e Vh{D)nBp, 

{H6) distiwo,M) <S = Cil + \\wop\\i,D + \\wh\\i,D)\\wo\\i,Dh^''+K 
where Bp = {w\w G V{D), < p}. 

Proof: Let Wi G Vh, = <p{wi) for i = 1,2, and <p = <pi — <p2- If setting w = Wi,W2 in (10.22) 
respectively and making subtraction, then we get 

Awh{4>,v) = bo{Wi - W2,Wi,v) +^boiW2,Wi - W2,v) - Awhi'^l - W2,V] 

-Awh{v,wi-W2), yv^Vh- 
Owing to Wh, is nonsingular and (10.13), the fohowing inequaUty is vaUd, 

2Q!o||0||i,z? < sup ^ < sup ii^^ii^ ^ 

< M(||wi||i,D + ||tU2||l,D + ||tf/t||l,D + l)||wi - W2||l,D- 

Furthermore by using the triangle inequality d < Hifo — Wh\\i d + ll^foHi d, we derive the 

(H5). 

Our task is now to prove (H6). We first know that 

dist(tuo,A^) = inf \\wo-w\\id 
weM 

< \\wo - (wop + cj){wop))\\l,D 
= \\wo - Wop - (f){W0p)\\l,D 
= l|l«Oq - 0(l«Op)||l,D- 

For any v G V{D), Equation (10.1) can be rewritten as, 

•Awh i'Woq, v) + Awh {wop, v) + bo{wo -Wh,Wo-Wh,v)- bo{wh, Wh,v) =< G^,v > . (10.24) 

Choosing w = Wop in (10.22) gives 

Awh{(t>{i"Op),v) = bo{wop,wop,v) - Awh{'t"Op,v) - Awh{v,wop)+ < Gn,v >, Vi; e T4(D). 

(10.25) 

Setting V e Vh{D) in (10.24) and Subtracting with (10.25). Taking into account (7.15), then we 
have 

Awhiwoq - 4>{WQp),v) = -bQ{wo - Wh,Wo -Wh,v) + bo{Wh,Wh,v) - bo{wop,WOp,v) 

= -bo{wo - Wh, Wo - Wh, v) + bo{wh - wop, Wh,v)+ bo{wop, Wh - wop, v). 
Since wo is a nonsingular solution of (7.1), and noting that (10.5) and (10.15), we get 

ao\\Woq - (t>[Wop)\\i^D < sup °- ii^jir— <I + II. 

vev(D) " '"'"^ 
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Our problem reduce to the estimation of /, //. It is easy to show that 

Aw AW o<,-4>iW Op), V) 
^ = \\V\Un 

_ Aw^iWog-4>(.Wop),Vop+Voq) 

~ ivopj:!), vm 

Aw AW a,- (t){W Op), Vo„) 

< sup ^r^^ 

Voqd Vh(D) 

_ -ha{Wo-Wh,Wu-Wu/V^ui) + hAWh-W^y,,:Wh,ra,,) + h,(Wp/W,,-Wop,Vo,) 

Voi^VhiD) 

1 3 

< - WhWo dW'^O - '^h\\lD + + ||Wfc||i,D)||l«Op - Wh\\l,D) 

< M(l + ||«;op||i,D + ||«'/i||i,d)I!«'o - whWq^dWwq - whWI^d-: 

jj _ Awo(Wog-(t>(Wop),V)-Aw^{Wog-(t>(Wop),V) 

~v!7iD) ii^ii^'- 

_ bo(Wo-Wh,Wog-4>{Wop),V)+bo{Wog-4>{Wop),Wo-Wh,V) 

~ P \\V\\l D 

vev(D) " 

< 2M\\Wo - Wh\\l,D\\woq - 0(t«Op)||l,Z3 

< 2MCh^\\wQq-(P{wQp)\\x,D. 

where in the fifth step estimation of I we adopt (10.16). Finally, noting that ( HI), we have 

<C/l2fe+S. 

So we prove (H6). □ 

Theorem 10.4. Assume that the assumptions {HI) ^ {Hi} for finite element space 14 are satisfied 
and Wh is the nonsingular solution of (10.3). Then variational problem (one step Newtonian 
iteration) 

J Findw^ e V{D)such that 

\ Awh{w*,v) =< Gn,v > +bo{wh,Wh,v), \/v e V{D), 
has a unique solution tu* and the following estimation are valid 



(10.26) 

\\w-w4i.D<Ch^''+', (10.27) 



where 

1, for B.C.I, 



^ ^ i, for B.C. 1 1 
Proof: Navier-Stokes equations (10.1) equals to 

Awh{wo,v) + bo{wo - Wh,wo - Wh,v) - bo{wh,Wh,v) =< G^,v >,Vi; e V{D). (10.28) 
Subtracting (10.28) from (10.26) leads to 

Awh{wo-w,,,v) + bo{wo-Wh,wo-Wh,v) =Q,\/v & V{D). (10.29) 
By applying (10.13) and lemma 10.1 we assert 

1 II 11/ AwAWa-W.^v) 

vev,AD) 

^ s^p — imr^ — 

^ -boiWo-Wh,Wo-W.,V) 
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For B.C.I., we have 

\bo{wQ - Wh,wo - w*,v)\ = I - bo{wo - Wh,v,wo - w^.)\ 

< M\\wo - Wh\\o.D\\Wo - Wh\\l,D 

Similarly, For B.C. II., 

\bo{wo - Wh, wo-w^,v)\< M\\v\\i,d\\wq - WhWlnWwo - Wh\\ln < MC\\v\\i,D\\wo\\l+i^Dh^''+^ ■ 
Thus the proof is completed. □ 
Remark 10.1. First, it is simple to show that 

Second, the variational problem (10.26) is still an infinite dim,ensional problem. We can apply the 
standard two-level finite element method on this problem(see Layton et al. [21-23] and references 
therein). 

Theorem 10.5. Suppose that the assumptions in theorem 10.4 satisfied. Vw is a finite element 

subspace with mesh parameter h* < h and satisfies assumptions (HI) ^ (-ff4) with integer m < k. 
If {wh*) is a Galerkin finite element approximation solution to (10.26), that is 



Find{wh*) G Vh*such that 

Awh{i"h*,v) = {Gr,,v) + bo{wh,Wh,v) Vi; e Vh* 



(10.30) 



Then the following error estimation holds 

\\W^ - Wh4hD < Ch^'^+'Wlw^lrn+l). (10.31) 

Proof:The proof is omitted. □ 

Combining theorem 10.4 and 10.5 leads to our final conclusion, 

Theorem 10.6. Suppose that the assumptions in theorem 10.4 o-nd 10.5 are satisfied. Then we 
have following estimation 

In particular, if choosing h* = then we have 

\\{W0-WH*)h,D<C{h'"'+'). 

Proof:The proof is omitted. □ 

Algorithm 1. Here we present the finite element approximation algorithm based on the approxi- 
mate inertial manifold, i.e., 

• Stepl: Solve the nonlinear problem (10.3) on the coarse grid with mesh size h, 

• Step2: Solve the linear problem (10.30) on the fine grid with mesh size h* . 

Remark 10.2. If we use the linear finite elem,ent method for (10.3) and (10.30), respectively, then 
on the two dimensional problem the following estimates hold 

\\{wo - Wh*)\\i,D < ch^ « c/^*^ 
where h* » hi . As we know that in Layton[19, Theorem 2], the result is 

\\{wo - Wh')\\i,D < ch"^ ~ ch*, 
with h* h'^. This shows that our results is much better than that in [19]. 
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A Appendix 



This section gathers most of the prehminary knowledge that will be required in this article. In 
subsection 1, wc focus on the the expressions of some physical and geometrical quantities in the 
new coordinates system; then in subsection 2, The Navier-Stokes Equation in the new coordinate 
system is derived. Finally, we consider the Gateaux derivative of the solutions of the Navier-Stokes 
equations with respect to the shape of blade. 



A.l Some Physical and Geometrical Quantities 

In order to simplicity, wc consider the 3D fluid flow in an flow passage in an impeller with rotating 
angular velocity us = (0, 0, w) around its axis, and the thickness of the blade is uniform. Let 
(x, y, z) be the cartesian coordinate system out of the impeller in the Euclidean space i?^, and three 
coordinate basis vector arc i,j,fc respectively. Furthermore, {r,9,z) be the cylindrical coordinate 
system attached to and fixed on the impeller, and er,eg,ez are three basis vector of this system 
respectively. Next we define a new coordinate system {x^,x'^,x^) through the following relations 

,1 = ^, x^=r, x^=^ = e-\0-e{x\x^)), 

■ = x2, z^x\ e = e^ + e{x\x^), ^ > 

where <d{x^,x^) be a smooth mapping from a bounded smooth enough subset D <Z B? into R, 
especially, (a;^, a;^, 6(a;^, x^)) denote an arbitrary point on the blade surface. Let e, be the basic 
vectors of this new coordinate system. The parameter ^ satisfies < ^ < 1 and obviously, ^ = const 
represent a surface 3j in B?, which can be obtained by a rotation of the blade through an angle 
of degree. 

Next, we present the following proposition. 

Proposition A.l. The covariant components Uap and contra-variant components a"^ of the met- 
ric tensor of the surface Ssj and the covariant components gij and the contra-variant components 
g^^ of the metric tensor of 3D Euclidean space are given by, respectively, 



aa/3 = 5a/3 + (x2)'e„e^, a = det{a„0) = 1 + {x^ (Qj + Q^^) , e„ = fl 

P" Va' Va' V" (A I 2) 

gap = dap, 53a = 5a3 = ,933 =e^?'^, 9 ^ det{gij) -^"^ 



where dj is the Kronecker symbol. Furthermore, we have the followings conclusions 
1. Rotating Angular Velocity w 

1 — j.i <,,2 _ n ,.i3 — 



CO = LO 



uP- = 0, = -w£-iei, 



2. Coriolis Forces 



= 0, = -2ruT\.{w, 6), (A./.4) 
C3 = 2c<;e-i(re2n(«),e) + 5^), 
3. Unite Normal vector to the Surface 3j 



n=-x''QJ^ae^ + (ea;2)-ii^e3, 
n« = -a;2e„/^/^, = {ex'^Y^^-^^ . 

4. Curvature Tensor of the Surface 9^ (Second Fundamental Form) 

611 = 7s(02(aii - 1) + x^Gii). 612 = 621 = ;^(eiai2 + x^Gia), 
bii = ^(©2(022 + 1) + x^Qii), b = det{bafi) = 611622 - 6?2> 



{A.I.5) 



{A.I.6) 
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where 



|ve|2 = e? + e»i. Ae = a«^e„^, Ae = 611 + 622, 6„^ = a„a^6, 

5. Mean Curvature and Gaussian Curvature of the Surface 

r K = b/a, 

\ 2H = ^[a;2(a226ii + 011622) - 2012612)62(2011022 + on - 022) - 26io?2> 
Proof: Firstly, from (A.I.I) we have 

X = x{x^,x'^, ^) = rcos6 = x^ cos(e^ + Q{x^,x^)) 
y = .T^,^) = r sin^ = a;^ sin(e^ + 6(a;\a;^)) 



Therefore, 



dx 

dx^ 
dy 

dz . 
3x1 - 



z{x^,z'^, ^) = x^ 



dx 



{A.I.7) 



{A.I.8) 



{A.I.9) 



— a;^sin^6i, 7— ^ = cos^ — a;^ sin^62, 

ox-' 



X cos ^61, 



dy 



dz dz _ Q 



dx^ 

' dx-^ dx^ 

where 9 = e£, + Q{x^,x'^). From (A.I. 10) wc get 

d{x,y,z) 



sin^ + a; cos ^62, 



dx 

F 

dy 
d^ 



■ —x^ sin^£, 
x^ cos 6e, 



(A7.10) 



d{x^,x'^,x^) 



= ex 



It is well known that 



Br = COS 6i + sin 6j, eg = — sin 6i + cos 9j, 
i = cos 9er — sin 6eg , j = sin 6er + cos 6eg , 



(AJAl) 



(y4./.12) 



Let 3?(.T-'^, .T^, ^) = x(x^,x'^,S,)i + y{x^ jx"^ + z(a:;^, a;^, denote the radial vector at the 
point P = (x^,a;^,^). Then the covariant basic vectors (e^.e.-j) and the contra-variant basic 
vectors (e",e'^) in the new curvilinear coordinate system {x/",^) arc given by, respectively, 



(A./.13) 



Ba = daU = daXi + d^yj + daZk, 63 = = if * + Ifj + if fc, 

ei = x^6iee + k = —x^ sindOii + x'^cos6Qij + k, 

62 = 92a;^ee + = (cos 6 — x^ sin 9Q2)i + (sin 6 + x^ cos 9Q2)j, 

63 — x^eeg = —ex"^ sin 0i + ex^ cos 9j, 

e' = (/*^ej, = e^j — e~"^6„e3, = — £"-^60,6^ + (r"e)~^oe3, 
e^=k, = cos 6*1 + sin6'_7, 

= — (re)^-^(sin6' + r62 cos6)i + {re)~^{cos9 — r62 sin0)j — e~^6ife. 

Inversely, we have 

Cr = 62 - £"^6263, e0 = {ex^y^e3, A; = ei - £-^6163, 
i = cos ^62 - (£~^cos^62 + (£a;2)"^ sin 6')e3, {A.I.IA) 
j = sin6e2 + ((£a;^)~^ cos6 — £"^62 sin^)e3, 

For any fixed ^, the mapping 

^{x\x^;0=x{x\x^,Oi + yix\x'^,^)j + z{x\x'^,^)k 

define a 2-dimensional surface SSj with single parameter ^ and the covariant components of metric 
tensor of is expressed as 

aafi = e^ep = Sc,^ + {x'^fOaQ^, a = det(o„^) = 1 + {x'^f{Ql + 6^) = 1 + |V6|2. 
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Meantime, the covariant components of the metric tensor of 3D EucHdean space in the coordi- 
nate system {x^jx"^,^) are given by 

Prom this and the expression (A.I. 13) it is easy to derive (A.I.I). 
Next, we consider the angular velocity vector. Obviously, 

which norm can calculated as 

= w • w = gijUj^u^ = gnu^u^ + g^s.Lo'^ lo'^ + 2gizu^oj^ 
= aii(w)2 + r2£2(w)2e-2@2 ^ 2£r2ei(w)2(-e-iei) 
= (a;)2(an + r^G? - 2r^el) = (w)^. 

Similarly, from the coordinate relation (A.I.I) and the definition of Sijk, the Coriolis force is 
formulated as 

C =2u: X w = 2eijk0o^w'^e^ = {2eijkUJ^w'^)g^'^em 

= 2s,,kt^^wHr/^e, + g'^e2 + g'^e^) = C^e^, 
Ci = 2{g'haki^w'' + g'^eak{-uje-^Qi)w'') 

= 2w(5"ie„i3ii;3 + ,g3ie3^2w2 - ff'^^eas^je^'eiu;^) 

= 2a;(0 - e-^Qiw^^ - ei32e~^Biw;2) = o, 

= 2{g^eak^w'' + g'^eak{-^e-^ei)w'') 

= 2uj{g''^eciaw^ +g^^e3i2w^ - e-^Qig^^eo^^pvjf^) 

= 2a;(£2i3W^ - e~^92w2e3i2 - e""^©ie23i^f 

= 2bj^{-w^ - £-192^2 - £-191^1) = -2ruj\l{w, 9), 

= 2{g'haktow'' + g'^eak{-uje-^Qi)w'') 

= 2(g"3^„i3wu;3 ^ ^33^g^2W^2 ^ g<^3eo,:if3{-uje-^Qi)w^) 

= 2a;(-£~i92£2i3M^^ + g^^£?.i2uP - e^^Qi{-e^^<d2£23iw^ - £~^^i£i32w'^)) 
= 2uj^{e-'^Q2W^ + g^^w^ + £-291(82^1 - Q^w'^)). 

Owing to the identity 

^33 = (r£)-2y,2 ^ ^-2(02 ^ 02)^2^ 

we have 

= 2rw£-^©2n(«;, 9) + 2uj{re)-'^v? . 

where 

n(u;,9) = £M;3 + 9aw". 

Therefore, 

2ijjxw = -2rwn(w;, 9)e2 + (2ra;£-i92n(u;, 9) + 2u){rs)-'^w^)e3] 
Finally, the contravariant components of Coriolis force in the new coordinate system is given by 
(71=0, C2 = -2wrn(u;,9), = 2w(£)-i(r92n(w.9) + ^f^). (A7.16) 
Next we consider the unite normal vector n of Q^. At first, it is well know that, 

n = r = ^(ei X 62) = n^t + Uyj + n^k = •n/Bi 

ei X 62 va 
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By virtue of (A.I. 13) and (A.I. 14), the above expression shows that 



i j k 

ei X 62 = (ei)a; (ei)j, (ei) 
(e2)x {e2)v (62) 

i j k 

—rQismO iQiCOsO 1 

cos 6 — r02 sin 9 sin ^ + rB2 cos 9 

= — (sin^ + r@2 cos 9)i + {cos 9 — r©2 sm9)j — r@ik 

From this we obtain the contra-variant components of n in the cartesian and the new coordinate 

system, 

rix = —^{sm9 + x^Q2 cos9), Uy = -^{cos9 — x^Q2 sin0), 
riz = -a;2©i/-v/a. 

= -x'^Qa/s/o., = {re)~^^, 

Now we calculate the curvature tensor of the surface Noting that 

1, , 1 

Oa^e = -^^{naeis + npea) = nsa/s = —j=^\ x ^2 • e^^, 

If the radial vector at the point P on 9^ is denoted by 3ft, then 

(923ft 



dx'^dx^ 



where Xa/3 = dad^x. Therefore 
\/aba0 



Xal3 Uap Zaji 

(ei),. {ei)y (ei)^ 
(62)0; {e2)v {e2)z 



Xal3 

— rGi sin 9 

cos 9 — r92 sint 



Vafi 

r0i cos 9 



sin 9 + rQ2 cos 9 



= - [{Xai3 sin 9 - ya/3 cos 6*) + r02 {Xai3 cos 9 + yai3 sin 6*)] 

= - [{Xa0 + r@2ya0) Sm9 + {rQ2Xa0 - Voip) cos9], 

Simply calculation from (A.I.8-A.I.10) shows that 

xii = = -x'^{&iism9 + elcos9), 

X\2 = —©1 sin 6' — x'^{Qi2 sm9 + 6162 cos0), 
a;22 = -262 sin 9-x'^ (822 sin 6* + 9^ cos 6*) , 

2/11 = a;2(0ii cos 6' - 9^ sm9), j/12 = 9i cos^ + a;2(6i2Cos^ - ©162 sm.9), 
2/22 = 2©2Cos6' + x2(©22Cos^- ©^ sin 6*), 

Zap = 0, 

Substituting (A.I.20) into (A.I. 19) leads to 

611 = 7;^ (02(aii - 1) + a;29n). 612 - 621 = ^{'diai2 + ^2912), 
h22 = ^(92(a22 + 1) + a;2922), b = det(6«/3) = ^>ii^>22 - 6^2, 
Finally, the mean curvature and the Gaussian curvature are calculated as 

I 2H = ^[x2(a22©ii+aii©22)-2ai2©i2)©2(2aiia22 + aii-a22)-2©iaf2, 
Those are (A.I.4)(A.I.5). 



{A.I.n) 



{A.I.IS) 



{A.I.IQ) 



{A.I.2Q) 



{A.I.21) 



{A.I.22) 



□ 



In the next place, we consider the Christoffel symbols and the covariant derivatives under the 
new coordinate system. 
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Proposition A. 2. Under the new curvilinear coordinate system {x",^), the Christoffel symbols 

and covariant derivatives are respectively given by 

= -rS2aOi:sOj, T^i^ = -erS2a'df3, 

rL=rL = r--'^2a+re2e„ T^s = -s^rS2a, rl3 = erQ2, 
and 

Vaw^' = l^-rS^Qalliw,e), 

^aw^ = s-Hx'')-^w^^a + e-^w^e^^ + iex^)-^a2an{w, e), 

= x^eS2an{w, e), v^w^ = 9^ + ^ + x^e2ii{w, e), 

Proof: Prom (A.I. 13) it follows that 

en = -x^{cos0Ql + sin6'eii)i + x^{-sm9Ql + cos6'eii)j, 
612 = 621 = (— sin^^Gi — x^(cos d<diQ2 + sin 0Qi2))i- 
+(cos6'ei +a;2(-sin6'eie2 + + cos 6'ei2))j, 

622 = (-2sin6'e2 - a;2(cos6'e2e2 + sin6'e22))i + (2 cos 6*6 1 + sin 616262 + +cos6l622))j, 
ei3 = 631 = —{re)Qi{cos6i + sinOj), 633 — —re^{cos9i + sin6'j), 

623 = 632 = — e(sin 6 + rQ2 cos 9)i + e(cos 6 — r62 sin 6)i, 

= e^en = fcen = 0, rf, = -x^Ol Tf, ^ ie)-\re20l + 0,,), 

r}2 = r^i = o, = r^i = -reie2, rl^ = rl,^ ire)-Heia22 + rei2), 

Tl^ = 0, - -rOl = (r£)-i[2e2 + r(e22 + rei)] - (r£)-i[62(l + 022) + re22], 

r}3 = rii = o, Tl,^Tl = -{rer'&u r?3 = ri^ = reie2, 
r^3 = r32 = o, ri3 = r^^ = -ree2, r^g = ri2 = r-ia22, 
r^3 = o, rl^ = -re\ ri3 = -r£62, 

This yields (A.I.23). The (A.I.24) can be obtain from ( A.I.23) and 

This ends the proof. □ 

A. 2 The Navier-Stokes Equation In the New Coordinate System 

Proposition A. 3. The Rotating Navier-stokes equations in the new coordinate system {x",^) can 

be written as 

' lf^ + ^ + ^=o, 

J\f"{w,p,e) := -i/Aw" + VaP + C"{w,u) - i^l"{w,e) + ^{-i^l'^{w,e) 

-e-^e^p) + NJ{w,p) +M^{w,p) = r, (A7/.1) 

AfHw,p, 6) := - e-ie„^ + C^w, - i^l^{w, 6) 

+ -§^{-vll{w,e) + (re)-^ap) + N^{w,p) + Afl{w,p) = f, 

where C{w,ij) is Coriolis forces defined in (A. 1. 4), and the other shortening symbols are definitely 
expressed as, respectively, 

' N^{w,w) =wf'^-r62all{w,eMw,e), Af^"{w,p)=w^^, 

< N^iw,w) =w'^^+s-'^w^w^ef3x + {re)-^Il{w,e){2w'^+r^e2li{w,e)), {A.II.2) 

^^^i{w,p) =w^^, 
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and 



l"{w,e) =-2reS2aQx^ + l^-^S2a + S2a{rAe-2ae2)n{w,Q) 
-[{rs)-\6cxe2 + 262aex) + e-^6„xAe]w^ - 2r-^52aw\ 

B^{e) =S2a[2{S2am^-rexeXa)], 



+ae2e2w^ + Bi{e)w'' , 

q{w,e) = ire)-^a^ -2s-^^Q0^ +2e-\r-^S2a - 
+e-i(re2|Ve|2- Ae)«;3, 
[ Blie) = (r£)-i[(r-i + raQ2e2)e, + 2Q2^), 



(A.IIA) 



e 



Ae = e„„ 



e 



611 + 622, |V6|2 = e2 + e2, 



Proof: Firstly, we derive the Trace Laplacian operator in the New Coordinate System, that is, 

{A.II.6) 



Aw°' ^ g^^ViVjw"" = Aw" + ^^^(ii;, 9) + 6), 
Aw3 = g^^ViVjW^ = Aw^ + d^lfiw, 6) + P{w, 6). 



In fact, we have 



Employing (A.I.23) and (A.I.24) we claim that 

= dp{^ - r<52a6^n(w;, 6)) + (r|^(5^ - r^p^S^„)Vxw- 

-e-i6;3(V3V^w« + V^V3W") + (r£)-2aV3V3w;" 



and 



= I [(re)-2a^ - {re)-H2al^{w, 6) - 2e-^Qp |^] + Qp52a gf? (rn(t«, 9)), 

Summing up the above results, we get 

Aw« = dpdpw°^ + J|[(r£)-2a^ - (re)-i,52an(u;,9) - 2e-^@f,'^] 
-52arA@n{w, 6) + hVxW + hVaw^ + hVxw^ + hVsw^ 

Aw'- + ^[(re)-2a%^ - {re)-'S2a.n{w, 9) - 2e-'Qp'^ + hw^ + hw^] 



+h 



[{re) 'a 

3 dx^ 



2„ dw" 



[-r52aQxh - reS2xl2 + {re)-^a2xh + r92/4]n(tz;, 9) 



+ [e-\Qxa - r-^Qx52a)h + r-H2M]w'', 



{A.II.7) 
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where some marked symbols are 

-{re)-'^a{-r£'^)52x5oca = r-^S2xSaa, 

h = -r$pSax + e-'ep{2Ty^x - t^a) + (^^)-'«(r^A - ^h^c^x) 

= -[{rs)-^{a2^Q0 + 62) + e-^AQ]Sc,x 

+e~^9^(2r~^a2/3^aA + r52aQxQi3) + {re)~'^a{-re52aQx - rs@2Sax) 

= -(re)-i(<5„Ae2 + ^2aeA)-£-i<5«AAe, 
73 = T^^-s-^QxT^s = -rsS2aex-e-^ex{-e^rS2a) = -2rsS2aQx, 
h = (re)-'ar^3 - s-^QpT^p = {re)-^a{-rs-'82a) - e-^Qp{-re52aQ0) = -r-^S2a, 

Therefore 

-rd2a<dxli - reSixh + {re)-'^a2xh + r@2h = {rAQ - 2aQ2)S2a, 
e-H&Xa - r-^Qx52a)h + r-H2ah = [r-2(2a - 3)62^ - 2rexexa]S2a- 

Substituting the above expression into (A.II.7), then 

Aw;" = Aw» + ^[{re)-^a^ - {re)-'52an{w,e) - 2e-'e^^ + hw^ + hw'] 

+r-i ^ - 2re52aQx ^ + [(r AO - 2ae2)52a]n(«;, 9) {A.ILS) 
+ [r-2(2a - 3)52a - 2rexQxa]haW'' . 



(A.II.8) can be rewritten in a splitting form, that is 

d /Q 



Aw" = Aw"" + -^Ifiw^Q) +l°'{w,e), {A.II.9) 



where 1°" (w , Q) , l'^ {w , Q) , {&) are formulated in (A.II.3). 

Our task is now to prove the second equality of (A. II. 6). Indeed, 

Aw^ = VfiVfivy^ - £~^<dp{V dV^w^ + VaV/sw^) + {rey'^aV^V^w'^ . 

An argument similar to the one used in proof of the first equality of (A. II. 6) shows that 

Aw^ = d0{Wf}W^)+Tl^Vpw^ + iTl0-T{^)Vpw^-T$^V3W^ 

+d^v^w^ + rl^Vpw^ + (rigj^A - ^^fs)Vxw^ - rl^Vaw^] 



and 



where (by (A.I.23)) 

>x-e-'e 

(re)-2ar^3 + e-^Qxi-rhS^x + 2rf J = 2r-^620 



■ Jl 


= i^x 




■ 1 


= ^1/3 ~ 


-rix- 




= ire)- 






— '-0C 


j + 2e- 



^ „ 



^0'- 0X 

■"3/3 = (^. 
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Some further calculations show that 

= dp{^ + {re)-\{5^2Qi3 + re^p)w" + a2pn{w, 6)]) 

-s-'e^dpi^ + ^ + re2U{w, e)) + f^[{re)-^aV3W^ - e-'Q^Vpw^] 

+ g|,((r£)-i[(5„2©/3 + re^p)w'' + a2pU{w, 6)] - s-'&fsC^ + re2li{w, &))) 
+e-^AQ{^+rQ2ll{w,Q)). 



Combining like terms, we get 

dfi{V0W^) - e-^QpdpV^w^ + ^[(re)-2aV3w3 - e-^QpVpw^] 

+ gf^(£-ie„0W«) + e-iAe(5^ + rQ2ll{w, 6)). 
On the other hand, by using (A. 1. 24), then 

+ [-r52xepJi + {re)-^a20J2 - re52xh + rOa J4]n(i«, 6) + (re)-i(2e2a - Ae<52a)t«" 

= Jidpw^ + J2di3w^ + d^iJ^w^ + JiW^) 

+ [£-V-2 + £-108262 - rs-^e2AQ]IL{w, 9) + (r£)-i(2e2a - Ae^aa)^^". 

Summing up the above conclusions, (A. 11.11) becomes 

Aw^ = Aw^ + ^[{rey^aVsw^ - e-'^QfiV^w^ - e'^^Qp^ + Jaw^ + JiW^] 
+ ^{e-'eapw'')+e~'A&{^ +j&2U{w,Q)) + Jidpw^ + J2djw^ 
+ [£- + £-106262 - r£-ie2A6]n(u;, 6) + (r£)-i(2e2c« - A652a)w"- 

(A//.12) 

Thanks to the expanded formula 

Jidpw^ + J2di,w^ = {re)-\620ex+rQpx)g^ + 

Hence, similarly, ^ 

Au;3 = Aw^ + ^lliw,e)+l^{w,Q), (A.//.13) 

where 6), 6) are expressed in (A. II. 4). This is the second expression of (A. II. 6). 
Our goal now is to consider the terms of the pressure. Actually, we have 

g'^^dfsp+g'^^d^p = 6''^dpp-e-^e^d^p, g^^^d^p+g^^d^p = -e-^QadaP+{re)-^ad^p. {A.II.U) 
By Using (A.I. 24), the nonlinear terms are formulated as 

w^Vjw" = w^Vpw" + w^Vsw" = w^^ +w^^ - rS2a^{w, e)U{w, 6) 
+(r£)-^in(t«,6)[2M;2 _^ [x^fQ2^{w,Q)], 
Combing (A.II.6) with (A.II.14)-(A.II.15) obtains (A.II.l). The proof is completed. □ 

A. 3 The equations for the Gateaux derivative of the solutions of NSEs 

In this section we consider the Gateaux derivatives of the solution of NSE with respective to the 
two dimensional manifold which is a portion of the solid boundary of the flow passage in impeller. 
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Proposition A. 4. Assume that there exists a Gateaux derivatives (w := ^^,p := of the 
solutions {w,p) of the Navier-Stokes equations (A. II. 13) with boundary conditions 

dw 

w\r, = 0, [ly— -pn]\r, = h, {A.III.l) 

where Tg and Ti are boundary dfl = U Fi, fl = D x [—1, 1]. Then {w,p) satisfy the following 
linearized Navier-Stokes equations with the corresponding homogenous boundary conditions, that 
is, 

-i^Aw" - vr{w, e) + C"(ii;, uj) + V„p - ^{lyl'^iw, 6) + e-^Qap) + N^{w, w) 

+N^{w,w) + N^{w,w) + N^{w,w) + R°'{w,p,e) = 0, {A.III.2) 
-uAw^ - ul^{w,Q) + C^{w,u)) - s-'^effdpp+ ^(-W^(tt),e) + {rs)-'^ap 
+N^{w, w) + N^{w, w) + N^{w, w) + N^{w, w) + B?{w,p, 6) = 0, 



and 



where 



w = 0, on r,n^ = a, _ (A III 3) 



R"{w,p)ii= -u%r)- -§^[^r, + s-^pr]o.] + ^T,-2rS^c.exIl{w,e)r]x, 

E?{w,p)r}= -v'^il+ ^[-'^■n+{re)-'^2r'^Qxprix] + ^'n-e-'^r]ccdcP [A.IIIA) 

+e~^w^w''r]xa + (rey^[w^{2w'^ + r'^e2^{'w,Q)) 

+n{w, e){d2xr^u{w, e) + r^e2w^)]m, 



{A.III.5) 



Proof:The Navier-Stokes equations (A.IFl) can be rewritten as 

^ + ^ + ^ ^0, 

J\f"{w,p,e)e^ + Af^{w,p,e)e3 = f^Ca+f^es. 
Set Gateaux derivative with respect witli 8 along any director 

rieW:=H^{D)n{H\D), with w\r, = 0, W^-pn]\r,=0} 

denoted by ^r]. Then from (A.IIL5) we assert 
^AA"(i^,p, e)e„J7 + ^^Af^iw,p, 6)6377 +AA"(«;,p, Q)^n + J\f^iw,p, 9)^7/ 

^^«(i«,p,e)e„ + ^^3(^^p^e)e3 + [Af-{w,p,Q) - f"]^ + [AfHw,P,e) - /3]^e3 = 0. 
Since Navier-Stokes equation (A.III.2), it yields 

f ^Af^iw,p,e)rj = ^U^(w,p,Q)v + ^U''{w,p,Q)wrj + ^Af''iw,p,Q)prj = 0, 
\ ^M^{w,p, e)»7 = ^M^{w,p, 6)77 + ^M^{w,p, e)wrj + g^"{w,p, 6)^77 = 0, 

{A.III.8) 

It is obvious that 

R"{w, 6)77 = ^N<^{w,p, 6)77, R\w, 6)77 = ^U\w,p, 6)77. 

□ 
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